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Abstract 



y—< The mechanism of magnetic flux maintenance on the solar surface in quiet regions on the 
O 

^ Sun is observationally investigated. There have been proposed two scenarios: One is by 
O 

Q the local interactions only in the surface level and the other is by the frequent replacement 
of the fluxes among different atmospheric altitudes. In order to understand the energy 

<N 

^ transport and budget through the different levels of the solar atmosphere, it is important 

^ to quantitatively study the surface magnetic processes since they are one of the major 

^ ^ drivers for the magnetic energy generation. The Solar Optical Telescope on board Hinode 
I 

O spacecraft makes, for the first time, such a study possible due to its stable and precise 

^ observation of the magnetic field. 

^ We investigate surface processes of magnetic patches, namely merging, splitting, emer- 

2 gence, and cancellation, by using an auto-detection technique. We find that merging and 

m 

'O splitting are locally predominant in the surface level, while the frequencies of the other 
(N 

T—t two are less by one or two orders of magnitude. The frequency dependences on fiux con- 
(N 

^ tent of surface processes are further investigated. Based on these observations, we discuss 



X 



a possible whole picture of the maintenance. Our conclusion is that the photospheric 



5_( 

magnetic field structure, especially its power-law nature, is maintained by the processes 
locally in the surface not by the interactions between different altitudes. We suggest a 
scenario of the flux maintenance as follows: The splitting and merging play a crucial role 
for the generation of the power-law distribution, not the emergence nor cancellation do. 
This power-law distribution results in another power-law one of the cancellation with an 
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idea of the random convective transport. The cancellation and emergence have a com- 
mon value for the power-law indices in their frequency distributions, which may suggest 
a "recycle of fluxes by submergence and re-emergence". 

We summarize the previous studies and the purpose of the thesis in Chapter 1. The 
description of data sets and method of the auto-detection code is explained in Chapter 2. 
The main results and discussion of the thesis are described in Chapter 3—5. 

In Chapter 3, statistical properties of detected patches are investigated. We use two 
data sets of magnetograms obtained by Hinode. One has higher time resolution but an 
intermediate duration (data set 1). The other has a long duration but lower time resolution 
(data set 2). The frequency distributions of flux content are power-law distributions with 
indices of — 1.79±0.18 (data set 1) and — 1.93±0.07 (data set 2). They are also consistent 



with the previous result by Parnell et al. (2009). The average lifetimes over the whole 



magnetic patches are found to be 17.3 minutes in data set 1 and 23.2 minutes in data 
set 2. These different values show that the obtained lifetimes depend on the temporal 
resolution of the data sets, i.e. they become shorter as measured with a better resolution. 
The lifetime of patches increases with the larger flux content and saturates at ~10^^ Mx to 
the value of 60 min. We investigate proper velocity of patches in data set 1. The averaged 
proper velocity is 1.2 km s~^. The flux dependence of proper velocity is found to be weak 
with a power-law index of —0.23. The frequency of splitting has a fiat distribution with a 
drop in the lower range of the flux content. It is interpreted that the actual frequency is 
independent of the flux content and that the obtained distribution is under the influence 
of the technical detection limit. 

In Chapter 4, we investigate frequencies of surface processes. It is found that oc- 
currence rates of emergence and cancellation are much less than those of merging and 
splitting. We found that probability distribution of merging has weak dependence on flux 
content, with a power-law index of 0.28. On the other hand, the frequency distribution 
of cancellation has a power-law distribution with an index of —2.48. 



Ill 



In Chapter 5, we present discussions of the results. First we summarize magneto- 



chemistry equation (Schrijver et aL, 1997), which describes number relationship between 
frequency distribution of flux content and those of magnetic processes. Splitting has a 
time-independent solution of a power-law distribution with an index of —2. Next, we also 
discuss the reason of steep slope of frequency distribution of cancellation. Our discussion is 
based on assumptions that the patch motions are driven by convection which has random 
flow direction with constant velocity and that the power-law distribution of flux content 
is predominantly maintained. The steep slope is naturally obtained there. The derived 



power-law index is very similar to that of the emergence (Thornton & Parnell, 2011). It 
suggests that large parts are re-emergences of submerged loops through cancellations. 
In Chapter 6, the conclusion of the thesis is given as follows: 

1) Frequency distribution of the flux content is maintained to a power-law distribution by 
merging and splitting on the solar surface. 

2) The frequency of cancellation can be interpreted as a result of collisions of patches 
under motions driven in random direction with constant velocities. 

3) Most of emergences are interpreted as re-emergences of submerged loops recognized as 
cancellations. 
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Chapter 1 



General Introduction 



Brief reviews of previous studies and the purpose of this study are summarized in this 
chapter. The main target of our thesis is magnetic field in quiet regions. It is affected by 
convective motion of plasmas on the solar surface. Magnetic structures on the solar surface 
are summarized in Section Magnetic processes i.e., merging, splitting, emergence. 



and cancellation, are reviewed in Section 1.2, The purpose and strategy of this thesis are 
mentioned in Section II. 31 



1.1 Magnetic Structures on the Solar Surface 

Magnetic field on the solar surface causes various energetic activities. Many researchers 



have been attracted in what properties it has and how it is structured (Figure 1.1). It is 
important not only as the energy source of solar activities but also as the only magnetic 
field on the stellar surface that we observe in the most detailed state. 

Magnetic field structure on the solar surface is classified into two regions, i.e. active 



regions and quiet regions (Figure 1.2). An active region is a strongly magnetized region 
that often contains sunspots and plages. Violent activities, solar flares, and ejections of 
plasmas (coronal mass ejections; CMEs) sometimes occur there. On the other hand, a 
quiet region was thought to be a moderate region without such activities. 
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However, recent high-resolution observations change our picture of quiet regions. One 
of the most noticeable results is that there is as large amount of magnetic flux content in 
quiet regions as in active regions. Recent observations by Hinode/X-Kaj Telescope (XRT) 
and SolarDynamicsObservatory (SDO)/ Atmospheric Image Assembly (AIA) show that 
there are a plenty of brightenings which is thought to be events releasing magnetic energy 
there. It is deduced here that a quiet region is not so moderate as we thought. There 
is another important character in quiet regions. The magnetic field is highly affected by 
convective motion on the solar surface, i.e. the plasma beta in quiet regions is higher 



than unity. Right panel in Figure 1.2 shows line-of-sight magnetic field in a quiet region. 
One can recognize that the field is organized as a collection of discrete units, here we call 
them magnetic patches. They are swept to the edges of convective cells. 



Magnetic field in a quiet region is divided into three categories, ephemeral region, 
network field, and internetwork field (same as intranetwork field in solar physics termi- 
nology). Our target is network field. Network field and internetwork field are ubiquitous 
structures in quiet regions. Network field is a magnetic field swept to the edge of super- 
granular cells. Convective motion is believed to have an important role in their formation 
and maintenance. The field strength of network field is from a few hundred G to ~ 1 kC 
Flux amount contained in each patch varies from 10^*^ Mx to 10^^ Mx (Martin, 1988, Wang 



et al. 1995). The lifetime varies from a few hours to a day (Thornton & Parnell 2011). 



The observed flux content of internetwork field varies from 10^^ Mx to 10^^ Mx (^Livingston 



& Harvey, 


1975 


Zirin , 


1985 


1987, 


Wang et al. 


1995 



field is in the middle of supergranular cells. The field strength of internetwork field is much 
weaker than that of network field. The smallest scale of internetwork field corresponds to 



the observational limit in the observation so far (Livingston & Harvey, 1975, Smithson 



1975 


Livi et al. , 


1985 


Martin , 


1988 


Wang 


1988 


Wang et al. 


1995) 



distribution of flux content is thought to continue below 10^^ Mx (Socas-Navarro et al. 



2004 , Manso Sainz et al. , 2004 , Khomenko et al. , 2005 , Dommguez Cerdena et al. , 2006 
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2011 February 13'' 8:00UT 201 1 February 13'^^ 8:00UT 
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Figure 1.1: Extreme Ultra Violet (EUV) image (left) and corresponding line-of-sight 
magnetogram (right) of the whole Sun. The EUV image is obtained by Atmospheric 
Image Assembly (AIA) on board Solar Dynamics Observatory (SDO). The magnetogram 
is obtained by Helioseismic and Magnetic Imager (HMI) on board SDO. The bright (hot) 
region coincides with the magnetized region. 
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Figure 1.2: Magnetograms of an active region (left) and a quiet region (right) obtained 
by Narrowband Filter Imager (NFI) on board Hinode. Left: Magnetogram of an active 
region. There is a negative sunspot in the middle of the field of view. Positive polarity has 
a less bundled structure than negative polarity. Right: Magnetogram of a quiet region. 
Both polarities are fragmented and seem to be affected by convective motion. There are 
some network structures (~ 20" — 40") in this image. 
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Rezaei et al. 


2007 


Sanchez Almeida 


2007 Harvey et al. 2007 


Orozco Suarez et al. 


2008 



Lites et al. , 2008). A quiet region looks a mixed-polarity region from network field and 



internetwork field. It is sometimes called magnetic carpet (Schrijver et al. , 1998, Parnell 



2001). The other one, ephemeral region, is a bi-polar region (Harvey & Martin, 1973 



Martin, 1988, Hagenaar, 2001), which is thought as a remnant of a small flux emergence 



Wang 


1988 


Harvey 


1993 



content than the other patches and larger proper velocity. The typical flux amount varies 



from (2 - 3) X 10^^ Mx to 10^° Mx (Harvey fc Martin 1973, Harvey et al. , 1975 Harvey 



1993, Chae et al. , 2001). The typical lifetime is from a few days to 4.4 hr (Harvey & 



Martin, 1973; Harvey et al. , 1975, Harvey, 1993). It significantly decreases with higher 



resolution. At the end of the lifetime, the association of the opposite polarities in each 
ephemeral region is lost due to the random transport by the convective motion. 



One question arises. What amount of magnetic flux is produced and transported? 
The first step to answer this question is to investigate total flux amount in quiet regions. 
It has been probed by the observations for a long time. However, we have not deter- 
mined the total flux amount because it increases as the spatial resolution of observation 
becomes better. One approach for this question is to investigate a frequency distribution 
of magnetic flux content in discrete magnetic patch, namely an investigation of scaling 
law. We study whether patches with smaller or larger flux amount are dominant in total 
flux amount from the shape of the distribution. Moreover, the distribution probably has 
the information for the flux maintenance mechanism. There are several studies on this 



issue (Wang et al. , 1995, Schrijver et al. , 1997, Hagenaar, 1999; Parnell et al. , 2009, Zhang 



et al. , 2010). But the distribution form is not settled for long time because it was differ- 



ent among above studies, such as an exponential distribution, a power-law distribution. 



and Weibull distribution. Parnell et al. (2009) gives one conclusion. It is a power-law 



distribution with an index of —1.85 ±0.14 between 2 x 10^'^ Mx and 10^^ Mx (Figure 1.3). 
There is a single power-law in the wide flux range spanning from large active regions to 
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small patches in quiet network. The collaborative usage of the whole Sun magnetograms 
obtained by Solar and Heliospheric Observatory (SOHO) /MDl and high-resolution mag- 
netogram obtained by Hinode/SOT enabled them to investigate in such a wide range of 
flux content. 

The total flux amount, $tot, is obtained as 

/'<)!'max 

$tot= / 0n(0)# (1.1) 

^ 0min 

where n{(j)) is a frequency distribution of flux content, (pmin and (pmax are the minimum 
and maximum flux content respectively. (Note that $ is called the 'total flux amount' 
but actually it is spatially averaged one over the area, namely its unit is Mx cm~^. In 
this thesis, we use this definition rule hereafter.) When n(0) = uq {(p / (f)o)~'^ , it becomes 



tot 



-7 + 2 



-7+2 / , \ -7+2 
^max \ / Ymin 



[1.2) 



In a case of 7 < 2, the frequency distribution is flat and patches with larger flux content 
is important for an estimation of the total flux amount. On the contrary, the frequency 
distribution steeper than 7 > 2 means that patches with smaller flux content is important 
for an estimation of the total flux amount. 

The active region has a typical flux range between 10^° - 10^^ Mx. The quiet region 
has a flux range below 10^° Mx. From n(0) = 1.0 x 10"^^ {(p/lO^^Mx)-^-^^ Mx^^ cm^^ 



from |Parnell et al.| ( |2009| , one obtains ^^J^ = 1.99 Mx cm'^ and = 1.50 Mx cm'^. 
The power-law index, —1.85, is not so different from —2, which cause an increase of total 
flux amount with higher resolution observations. 

Next question arises at this point. How the frequency distribution is achieved and 



sustained? Parnell et al. (2009) suggest two scenarios. One is that frequency distribution 
on the solar surface represents distribution of flux tubes generated by dynamo mecha- 
nisms in the convective layer. The other is that it is achieved and sustained by surface 




Figure 1.3: Frequency distribution of flux content reported by Parnell et al.(2009). 
Three different data sets are used for this plot. They are different in solar cycle-phase and 
dynamic range of observation. The obtained distribution is a single power-law distribution 
with an index of —1.85 ± 0.14 from large active regions to granular-scale patches. 



7 



1.2. SURFACE PROCESSES ON THE SOLAR SURFACE 



activities, namely emergence, merging, splitting, and cancellation of magnetic patches on 
the photosphere. The key investigation to distinguish them is to quantify the effect and 
amount of these magnetic activities. 



1.2 Surface Processes of Magnetic Patches on the Solar 
Surface 



The magnetic field near the surface is transported and entangled by the convective motion 
of the plasma. The appearance of such modulation is observed as various activities in the 
line-of-sight magnetic component. There are four surface magnetic processes to change a 



frequency distribution, namely emergence, merging, splitting, and cancellation (Schrijver 



et al. , 


1997; 


Parnell 


2001 



activities between two patches. Note that there are diffusions and appearances of the 



unipolar patch in observational data, which do not correspond any of them. Lamb et al 



(2010) reports that an unipolar appearance in intermediate resolution data corresponds to 
merging of small patches in high resolution data. It deduces that these unipolar process 
may be caused by the surface processes involving undetected patches. 

Flux emergence is a divergence of opposite polarities in line-of-sight magnetogram. 



The diverging velocity of polarities are reported by some authors (Brants, 1985, Zwaan 



1985 


Otsuji et al. , 


2007 


Otsuji et al. , 


2010 



early phase and decreases to ~ 1 km s ^. It is much larger than that of convective motion 
on the solar surface, ~ 1 km s^^. The emergence contributes to the increase in total flux 



amount. It is interpreted as a flux ascension from below the photosphere (Zirin, 1972 



Zwaan, 1987). Figure 1.4 (c) shows a schematic picture of it. Downfiow is observed at 



the foot point by Doppler velocity observation (Zirin, 1972). It is interpreted as dropping 



material by the gravity of the Sun. Upfiow is also observed in the middle of the separating 



polarities (Bruzek, 1969; Zirin, 1972). These signatures support a picture of ascending 
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(c) Emergence 




(d) Cancel I at ion 




Figure 1.4: Schematic pictures of magnetic activities between two patches (left column) 
and corresponding 3D pictures (right column), (a) Merging. Two patches converge and co- 
alesce into one patch, (b) Splitting. One patch is divided into two patches, (c) Emergence. 
Two opposite polarities with equal amount of flux content appear, (d) Cancellation. Two 
opposite polarities converge with each other. One of them disappears. The other loses its 
flux content or disappears. 



9 



1.2. SURFACE PROCESSES ON THE SOLAR SURFACE 



flux tube from below tlie pliotospliere. 



Flux cancellation is a convergence and a disappearance of magnetic fluxes of positive 
and negative polarities in line-of-sight magnetograms. This phenomena is observed in both 



active regions (Zwaan, 1978, Martin et al. , 1985, Zirin, 1985) and quiet regions (Harvey 



fc Harveyj [T976| |Livi et al.j [T985| [Wangj pMj [Harveyj p:985l [Harvey fc Z^i^ p^93l 



Harvey, 1996, Schrijver et al. , 1998). Two physical models, U-loop emergence and fi-loop 



submergence, are proposed (Zwaan, 1987). Figure 1.4| (d) shows the fi-loop submergence 
model. Harvey et al. (1999) investigate timing difference of cancellation in different two 
layers: One should observe cancellation earlier in the lower layer than in the upper layer 
in U-loop emergence and vice versa in i7-loop submergence. They obtained that most 
of cancellations (20 of 45 cancellations, 44%) are interpreted as an f2-loop submergence. 



Some recent studies report the Doppler velocity in canceling region ( Yurchyshyn & Wang 



2001 


Chae et al. , 


2002 


Kubo et al. , 


2010 


lida et al. 


2010; 


Chae et al. 


2010 


)■ 


Chae et al. 



(2002) report downflow at the horizontal field connecting canceling patches near active 



region. 



lida et al. (2010) report time evolution of downflow around cancellation site in 
network field. The lifetime of downflow is significantly larger than velocity fluctuation in 
surroundings. These papers report same downward Doppler velocity, 1 km s^^. As for 
more tiny cancellation in granular scale, significant signature of downflow has not been 



found so far (Kubo et al. , 2010). The other supporting material for i7-loop submergence 



is that cancellation in a quiet region often coincides an overlying bright point in EUV 



and X-ray image (Harvey, 1996). This suggests that the photospheric cancellation is a 
manifestation of the reconnection and corresponds a fieldline interaction in the upper 
atmosphere (Priest et al. , 1994, Parnell et al. , 1994a|[b Schrijver et al. , 1998; Litvinenko 



1999). Some papers report result of numerical calculation including convective effect 



(Cheung et al. , 2010). They found that cancellation takes place in a numerical box 



and it is an r2-loop submergence. On the other hand, Zhang et al. (2009) report the 



blueshift motion around the cancellation sites from results of the spectropolarimetric 
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observation. There are no statistical investigations of the direct measurement of Doppler 
velocity around cancellation sites so far. It is still an open question whether cancellation 
is an fi-loop submergence or a U-loop emergence. 

Emergence and cancellation have direct relation with flux product and flux loss in line- 
of-sight magnetogram. The flux replacement time scales by these processes are defined 
as 



''"emrg "^t 



9$t 



tot 



dt 



emrg i 



dt 



-1 



:i.3) 



:i.4) 



where, d^tot/dt is flux change rate of each activity. The previous papers show they 
drastically decrease from several days to several hours as spatial resolution becomes higher 



(Martin et al. , 1985, Schrijver et al. , 1998, Hagenaar, 2001). It is not clear why they 



decrease with higher resolution. 



Figure 1.4 (a) and (b) show schematic pictures of merging and splitting. They do 
not change total flux amount on the solar surface but change the flux content of patches. 
There are few reports for merging and splitting compared to those of emergence and 
cancellation. 

Some analytical models are suggested. Magneto-chemistry (M-C) equation is useful in 



this point. It is suggested by Schrijver et al. (1997). It describes a relationship between 



frequency distributions of surface processes and that of flux content. It is written as 



dt 



"OO 



H — / n±{x)n±{(f) — x)l±{x, (J) — x)dx — n±{(f)) / n±{x)l±{(l), x)dx 
2 Jo Jo 

/•oo rcf) 

+ 2 / n±{x)k±{(j), X — (f))dx — n±{(f)) / k±{x, (j) — x)dx 



'0 

oo 



+ / n±{(j) + x)nz^{x)m±{(f) + x,x)dx — n±{(f)) / nzf{x)m±{(j),x)dx 
'0 Jo 



[1.5) 
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where n(0) is a frequency distribution of flux content, S{(j)), l{x, y), k{x, y), and m(x, y) are 
functions representing probability density distributions of emergence, merging, splitting. 



and cancellation respectively (See Section 5.1 for more detailed explanation). There are 



five unknown functions in this equation. What we have obtained from observation is 
n{(j)) for time independent case. Some authors found sets of [S{(l)),l{x,y),k{x,y),m{x,y)] 



to make a time-independent solution of n(0) with some assumptions. Schrijver et al 



(1997) found one particular solution for an exponential distribution. They assumed 1) 



re-appearance of submerged flux by cancellation, 2) /(x, |/)=m(x, ?/)=constant, and 3) 



fc(x, y)=constant. another set was found by Parnell (2002). They assume re-appearance 



of submerged flux through cancellation. Probability density distributions are written as 

k 

k±{x,y) = ^f±{x)f±{y)g±{x,y) (1.6) 
3k 

l±ix, y) = —f±ix + y)g±{x, y) (1.7) 
n± 

I ^ fco f±{x-y)f±{y) , . , . 

m±{x,y) = — 9±{x-y,y) (1.8) 

n± Uiy) 

where k^ is a constant, n±{(j)) = n±/A f±{(j)), g{x,y) is an arbitrary function. Moreover, 
they assume 

/ N exp{-{x + y)/r]±) 

9±{x,y) = r ( I ( — , — 
f±{x)f±[y)f±[x + y) 

The deduced frequency distribution of emergence and cancellation becomes exponential 
in this solution. These solutions satisfy the frequency distribution of flux content in their 
study. 

We have not found a solution for a power-law distribution demanded from the obser- 



vational result by Parnell et al. (2009). It is difficult to obtain a proper solution because 
we need at least three assumptions to obtain all frequency distributions even in a time- 
independent solution. One approach is to measure the frequency distribution of processes 
directly from observation. The difficulty of this approach is that it needs a statistical in- 
vestigation, namely that we need to look huge amount of data or develop an auto-detection 
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code for surface processes of magnetic patches. Frequency distribution of emerging flux 



has been investigated by several authors (Harvey et al. , 1975, Harvey & Zwaan, 1993 



Harvey, 1993 Title 2000 Chae et al. 2001; Hagenaar, 2001 Hagenaar et al. , 2003 2008 



Thornton & Parnell, 2011). Thornton & Parnell (2011) gives one conclusion. They found 



a power-law distribution with an index of —2.69 spanning from active region to inter- 



network field by using Hinode/SOT magnetogram and the previous result (Figure 1.3) 
There are no reports of frequency distributions of other processes. 



1.3 Purpose of the Thesis 

A quiet region is thought to have an important role in flux and energy transport on 
the solar surface because it has a plenty of flux budget and a large number of energetic 



activities. Parnell et al. (2009) reported that frequency distribution of flux content has a 



power-law form with an index of —1.85 ± 0.14. 

Our main purpose is to investigate what makes a power-law distribution. The strat- 
egy of this study is that we investigate the frequencies of magnetic activities, namely 
emergence, splitting, merging, and cancellation, based on observations and construct a 
model from the result. The auto-detection code is needed and is developed in this study 
to investigate frequencies of surface processes of magnetic patches. Our target is network 



field. One reason is that most flux in quiet regions is contained in network field (Martin 



1990, Schrijver et al. , 1998). The other reason is that the lower limit of statistical study 



with recent high-resolution data, namely magnetograms obtained by Hinode spacecraft, 
is ~ 10^^ Mx that is good enough for network field while is insufficient for internetwork 
field. 

The thesis has 6 chapters. The detail of the data set and method of analysis are 
explained in Chapter 2. The description of data set and our auto-detection method 
(namely our definition of patches and activities) are explained there. Detail description 
of M-C equation is summarized in Chapter 3. The forms between frequencies of surface 
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\ 

_l I I I I I 

10"^ 10^ 10^ 10^ 

Flux (x10^' Mx) 

Figure 1.5: Frequency distribution of flux emergence reported by Thornton & Parnell 
(2011). They plot a frequency distribution of emergence on flux content obtained by their 
detection-code (indicated by B.C. and T.C in the figure) and previous papers (indicated 
by the first character of authors and published year) together. Their fitting results in a 
single power-law distribution with an index of —2.69. 
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processes of magnetic patches and physical quantities are deduced there. We introduce a 
detection limit for statistical study, above which we observe all phenomena in statistical 
sense. Flux dependence of various properties of magnetic patches, size, flux content, 
proper velocity, and lifetime, are investigated in Chapter 4. Frequencies of magnetic 
activities are investigated in Chapter 5. We suggest one qualitative model in the discussion 
there from the observational results. The summary of the thesis and future works are 
summarized in Chapter 6. 
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Chapter 2 



Data Description and Method of 
Analysis 



We show data description used in the thesis and our method of detection of patches and 
surface processes of them in this chapter. Two sets of magnetograms obtained by Hinode 
spacecraft are used in this study. The instrumentation of Hinode spacecraft and data 



description are shown in Section 2.1 Our method of detection is explained in Section 2.2 



2.1 Data Description 

2.1.1 Nal Di Magnetogram Obtained by Hinode /NFI 

Two sets of hne-of-sight magnetograms obtained by the Narrowband Filter Imager (NFI) 
of the Solar Optical Telescope (SOT) on board Hinode spacecraft are used in this study. 
The reason why we use Hinode/NFl magnetograms in this study is that the observed range 
of flux content is wide enough for an analysis of network field, which is not accomplished 
with the other instruments. For example, the detection limit of MDI high-resolution data 



is ~ 10^* Mx (see Figure 1.3). We need wider range of flux content for a investigation 
of scale dependence of the network field (< 10^^ Mx). We summarize the description of 
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magnetograms obtained by Hinode/NFl in this section. 



Hinode is a Japanese mission of spacecraft observing the Sun (Kosugi et al. , 2007). It 



is developed by ISAS/JAXA, with a collaboration of NAOJ, NASA, and STFC. Figure 2.1 



shows Hinode spacecraft. There are three instruments on board, Solar Optical Telescope 



(SOT; Tsuneta et al. 2008 Ichimoto et al. 2008 Shimizu et al. 2008; Suematsu et al. 



2008), Extreme Ultra Violet Imaging Spectrometer (EIS; Culhane et al. , 2007), and X- 



Ray Telescope (XRT; Golub et al. , 2007). It is on a circular orbit with an altitude of 
~ 680 km. Hinode can continuously observe the Sun nine months a year. 



SOT is the largest telescope which has a Gregorian optics with 50-cm aperture. Its 
diffraction limit is 0.2" — 0.3" for the range of 3880-6302 A. The solar ray is divided 
to four paths in SOT, namely NFI, the Broadband Filter Imager (BFI), the Spectro 
Polarimeter (SP), and the Correlation Tracker (CT). NFI, which is used in this study, can 
record the filtergram image, dopplergram, and polarization image of some photospheric 
and chromospheric hues selected from Mgl (5172. 7A), Fel (5250. 2A, 5247. lA, 5250. 6A, 
5576.1A, 6301.5A, 6302.5A), HI (6562.8A), and Nal (5895. 9A). The full field of view of 
CCD camera is 328" x 164". The pixel size of CCD camera is 0.08". The actual observation 
type is set to meet the demand and follow the restriction of data telemetry. 

We use the circular polarization images of Nal Di resonance line, whose center is 
5896A, in this study. The formation height of this line is thought to be from the upper 
photosphere to the lower chromosphere. The Lande factor of this line is 1.33. The circular 
polarization images shifted from the line center are taken by SOT/NFI in our data set. 
They are converted to the magnetograms by considering the Zeeman effect. We employ 
the weak field approximation in this conversion. The detail of it is explained in the next 
chapter. The typical detection limit of longitudinal component of magnetic field strength 
is ~ 11 G depending on signal-to- noise ratio. 
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Figure 2.1: Outlook of Hinode spacecraft from Tsuneta et al. (^2008|. 



2.1.2 Description of Data Sets 



We use two data sets of circular polarization images of Nal Di resonance line. One has 
a high time cadence but the observing period is relatively short. The other has a long 
period but the time cadence of it is relatively low. From here, we call the former data set 
as data set 1 and the latter as data set 2 (Table 2.1). 

The time period of data set 1 is from 0:33UT to 4:08UT on 2009 November ll*'^. The 
time duration is 3 hour 39 minutes. Hinode observes a quiet region near the disk center. 
It obtains intensity signal (Stokes-I) and circular polarization signal (Stokes-V) at the 
two wavelength points shifted from the center of Nal Di resonance line by ±160mA. The 
magnetogram and dopplergram, which are calculated from these signals on the spacecraft, 
are downlinked to the ground. The center of the SOT view moves from (—14.1", 19.0") 
to (20.2", 19.4") in heliocentric coordinates during this period. The data are summed by 
2-pixels in both x and y directions. The pixel size of images becomes 0.16". Field of view 
in this data set is 113" x 113". Some network cells, which have a typical size of 20" — 40", 



are included in the field of view (Figure 2.2). The total number of images in this data set 
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Figure 2.2: Example of Nal Di magnetograms in data set 1 after the pre-process explained 



in section [2.2.2[ This magnetogram is taken at 2:21UT on 2009 November 11, which is 
the middle of the observational period. 

is 215 during the whole observational period. There are sometimes data missing lines in 
some images, which are lost during data transferring. We remove such images to obtain 
homogeneous images. The total number of polarization images is 199 after the removal. 
Time interval between consecuting magnetograms is ~ 1 minute. It sometimes becomes 
2 minutes due to our removal of data. 

The time period of data set 2 is from 10:24UT on 2008 December 30 to 5:37UT on 2009 
January 5^^. The time duration is 115 hour 13 minutes. Hinode tracks a quiet region near 
on the solar equator during this period. SOT/NFI obtains circular polarization signal at 
only one wavelength point shifted from Nal Di resonance line by 140mA. The circular 
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Data Set 1 


Data Set 2 


Observing Line 


Nal Di 5896 A 


Nal Di 5896A 


Pixel Size 


0.16" (2x2 summed) 


0.16" (2x2 summed) 


Observational Period 








2nnQ Nov 11 n-sniiT 


2008 Dpp 30 10-24TIT 


(F,nr\] 


2nnq Nov ii 4-nqTiT 


2009 Tan 5 5-37TIT 








I o L ai h j 






\£jl]\1 j 




^^uyo.y , — u.y j 


Duration 


3hr 39min 


115hr 13min 


Number of Images 


199 


1642 


Time Cadence 


~ 1 minute 


~ 5 minutes 


FoV 


113" X 113" 


94.9" X 91.4" 


Obsevation Type 


Stokes IV+DG 


Stokes-V/I 


Wavelength offset 


+/-160mA 


+140mA 



Table 2.1: Summary of data set description. 

polarization degree is calculated on the spacecraft. Note that the evaluation of magnetic 
flux density from this cirular polarization degree contains an uncertainty with an influence 
of the Doppler shift which is not measured in this data set. The center of SOT view moves 
from (—520.9", —9.3") to (695.9", —6.9") in solar coordinate during this period. The data 
are summed by 2 pixels in both x and y directions. Then the pixel size becomes 0.16". 
Full field of view in this data set is 113" x 113". The studied field of view is reduced to 
94.9" X 91.4" due to the removal of image edges for the correction of spacecraft jittering. 



As same as data set 1, some network cells are included in the field of view (Figure 2.3). 
Total number of images in the data set is 1642 after removing the data which have missing 
lines. Time interval between consecuting magnetograms is ~ 5 minutes and sometimes 
~ 10 minutes. 



2.2 Auto-Detection Algorithm 

We define four surface processes of magnetic patches and explain our method for identifi- 



cations of them in this chapter. Figure 2.4 represents schematic pictures of it. The basic 
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Figure 2.3: Example of Nal Di magnetograms in data set 2 after the pre-process explained 
in the section 2.2.2 This magnetogram is taken at 8:01UT on 2009 January 2, which is 
the middle of observational period. 
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Figure 2.4: Schematic picture of our patch tracking and detection f surface processes. 
(1) Detection of patches with a clumping algorithm. (2) Track of patches by examining 
overlaps. (3) Detection of merging and splitting by examining overlaps. (4) Detection of 
cancellation and emergence by making pairs of flux change events. 

chart is as follows: 1. Detection of patches with a clumping algorithm. 2. Track of patches 
by examining overlaps. 3. Detection of merging and splitting by examining overlaps. 4. 
Detection of cancellation and emergence by making pairs of flux change events. 

2.2.1 Preprocessing of Data Sets 

We have to convert polarization signal to magnetic field strength in data set 1. The 
dark current and the flat field of the CCD camera are calibrated by using the proce- 
dure fg_prep.pro included in the SolarSoftWare (SSW) package. We apply the Zeeman 
weak-field approximation to obtain the line-of-sight magnetic field strength from the cir- 
cular polarization (CP hereafter) signal and it is thought to be applicable to most of 
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network field (Landi DeglTnnocenti & Landi DeglTnnocenti 1973 Landi DegFInnocenti 



& Landolfi, 2004; Asensio Ramos, 2011). The proportional relationship stands up in this 
approximation. The condition for its approximation is that the effective wave shift of 
circular polarization by Zeeman effect is smaller that line broadening effect, which can be 
deduced as 



iTimc l2kBT 2 



where m is the electron mass, c is the speed of light, g is the effective Lande factor of the 
line, Aq is the wavelength of the line center, e is the electron charge, A;^ is the Boltzmann 
constant, T is temperature, M is the mass of Na atoms, and v^^ic is the microturbulence. 
We substitute Aq = 5896 A, ^ = 1.33, T = 6000 K, and timic = 1 km s~^ for this case 
and obtain the condition, B< 1800 G. This is satisfied in quiet regions. The conversion 



coefficient is determined by comparing the SOHO/MDl magnetogram (Scherrer et al. 



1995) and circular polarization in one image. We spatially smear the SOT data to two 
times the MDI pixel size (Dr. R. A. Shine, private communication) and make a linear 
fitting between CP in the SOT and the magnetic field in the MDI. Because the pixels 
in the SOT data, which are out of linear range, are too weak or too strong , we make 



the fitting in the range from 30 G to 100 G as a magnetic field strength. Figure 2.5 and 



Figure |2.6| shows scatter plots of signals obtained by SOT data numbers and magnetic 
field strength obtained by the MDI in data set 1 and 2. We obtain 9067.98 G DN'^ and 
1.83191 G DN~^ as a conversion coefficient respectively. 

Some corrections are done before the detection. Due to the change of the viewing 
angle into the solar surface in the long observing duration, line-of-sight magnetic field 
does not correspond to magnetic field perpendicular to the solar surface. We project the 
observed line-of-sight magnetic field to perpendicular field to the solar surface assuming 
that the magnetic field is perpendicular to the solar surface. There are sudden changes 
of field of view when the position of tip-tilt mirror hits the stroke limit due to the long 
observational period of data set 2. We remove such change of field of view by making a 
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SOT CP vs MDI Magnetic Field 




-0.02 -0.01 0.00 0.01 0.02 

SOT CP 

Figure 2.5: Scatter plot of CP obtained by Hinode/SOT and the magnetic flux obtained 
by S OHO /MDI in data set 1. Asterisks show the averaged SOT CP signal corresponding 
to 1 pixel of MDI. The horizontal bars indicate the minimum and maximum value of the 
SOT CP signal in each MDI pixel. The solid line shows a result of a linear fltting, whose 
slope is 9067.38 G DN~^. The fitting range is from 30 G to 100 G in absolute magnetic 
field strength. 
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SOT DN vs MDI Magnetic Field 
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Figure 2.6: Same plot as Figure [275] but for data set 2. The slope of fitting line is 1.83191 
G DN-i. 
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correlation of magnetic field between consecutive images. Then both data sets are rotated 
to the position at certain time, 2:03UT on November 11*'^ for data set 1 and 20:35UT on 
January 1** for data set 2, when the regions are near the disk center. There is sometimes 



a column wise offset of CP in NFI (Lamb et al. , 2010), whose reason is not known. We 
remove it by subtracting an average value of each line from all pixels in the same line. 
Magnetograms are averaged over 3 continuous images and 3 pixels for smoothing. Without 
these smoothings, the auto-detection does not give the stable results owing to the flipping 
noise around the threshold level. Note that these smoothings deduce the time-resolution 
of data set, namely 3 minutes for data set 1 and 15 minutes for data set 2. 



2.2.2 Detection and Tracking of Magnetic Patches 



Magnetic patches are detected in this step. We adopt nearly same definition as Parnell 



(2002) in this study. We use a clumping method for the detection of patches: each 
patch is picked up as a clump of marked pixels having magnetic strength beyond a given 
threshold. We set a signal threshold in this method. It is obtained by fitting the histogram 
of the signed magnetic field strength per pixel by Gaussian function in each magnetogram. 



Figure [27f] shows an example of histogram and fitting result in one image of data set 1. The 
horizontal axis indicates magnetic field strength and the vertical axis indicates number of 
patches per bin. Histogram represents the observational result and dotted curve represents 
the fitting results. Two vertical red lines indicate 2(T, which we use as the threshold. The 
actual value of a is ~ 5 G pixel ~ 6.8 x 10^^ Mx for each image. This value is close to 



that of Parnell et al. (2009), who report 7.0 x 10^^ Mx. We set a size threshold to pick up 
network field in addition to signal threshold. One scale on the network field is granular 
size. It has a typical scale of 10^ km (~ 9 pixel size in SOT/NFI). We pick up magnetic 
patches with sizes beyond 81 pixels for focusing our analysis on the network magnetic 
field. The validity of this choice is demonstrated in Figure 2.8[ which shows two-leveled 



magnetograms with different size threshold. As shown in this figure, one can pick up only 
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Figure 2.7: Histogram of the magnetic flux density in one magnetogram of data set 1. 
The histogram represents the observational result and the red dotted curve represents 
fitting result. Two vertical red lines indicate la. 



the network fields without internetwork fields by adopting 81-pixels as a threshold. 

The reasons why we adopt the clumping method are that it is simple and many of 
previous studies employ it. One of the limitations of this method is that we cannot 



investigate the inner structures of the patches. Note that DeForest et al. (2007) reportd 



that there are differences of the obtained frequency distribution of flux content among the 
detection method below the range of 1.5 x 10^*^ Mx with S'OifO/MDI data. We expect 
that there is no difference in the analysis range (> 3.1 x 10^'') because the resolution of 
Hinode/SOT is better than that of SOHO/MDl by more than one-order of magnitude in 
the flux range 

We track the motion of magnetic patches in consecutive images after the detection of 
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Figure 2.8: Two-leveled magnetograms taken at 2:21 on 2009 November 11*" with dif- 
ferent size thresholds. Size thresholds are set as (a) 9 pixels, (b) 36 pixels, (c) 81 pixels, 
and (d) 144 pixels. 



CHAPTER 2. DATA DESCRIPTION AND METHOD OF ANALYSIS 28 



them. Patches are marked as identical when they have a spatial overlap in continuous 



images (Hagenaar, 1999). The typical travel distance of magnetic patches in the data 
interval (~ 1 minute) is up to nearly 1 pixel size. However, a typical travel distance of 
patches becomes more than 1 pixel with a time interval more than 1 minute. We set 
margin to avoid this problem in data set 2. The margin is set considering the typical 



velocity, which we set 2 km s~^ from twice of a typical patch motion (See Figure 3.5), and 
pixel size. It becomes 5 pixels. More than one patch in a previous image often have spatial 
overlaps with one patch in a consecutive image and vice versa in a recent high-resolution 
magnetogram. To clear up this problem, we set two conditions when tracking patches. 
First, we check spatial overlaps from a patch with larger flux content. This is based on 
the concept that a smaller patch has a greater tendency to fall below the detection limit 
of the analysis by splitting and cancellation. Second we select a patch with the most 
proximate flux content in case of overlaps of more than one patch. Tracking paths of 



detected patches become unique with these conditions. Figure 2.9 show a tracking paths 
of center of detected patches in data set 1. 



2.2.3 Merging and Splitting 

Merging is a process where more than one patch of the same polarity converge and coalesce 
to one patch. A merging event is defined as a feature that satisfies the following conditions: 
(1) that there are two or more parent patches in a previous magnetogram overlap one 
daughter patch in the consecutive magnetogram, and (2) that more than one of the 
concerned patches in the previous magnetogram disappear in the time interval. Figure 



2.10 a) shows an example of the mergings. There are two separated positive patches at 
2:49UT. They gradually converge and coalesce into one patch at 3:02UT. They stayed to 
be one massif even at 3:08UT. 

Splitting is a process where a single patch is divided into more than one patch. A 
splitting event is defined as a feature that satisfies the two following conditions :(1) that 
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Tracking Result of Data Set 1 
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Figure 2.9: Example of tracking result in data set 1. The observational time is 0:30UT 
on 2009 November 11 The red (blue) circles indicate positive (negative) patches detected 
in this image. The black (white) lines show paths of center of detected patches. 
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there are one or more daughter patches in a magnetogram overlapping one parent patch 
in the previous magnetogram, and (2) that more than one of the concerned patches in 
the latter magnetogram appear in the time interval. Figure 2.10[ b) shows an example 
of the splittings. There is one negative patch at 1:50UT. A dip shape of the outline is 
formed between north-east and sourth-west parts of the patch at 2:06UT. It grows at 1:55- 
— 2:07UT and the patch finally splits into two negative patches at 2:10UT. The distance 
between them continues to be larger at 2:10— 2:14UT. 



2.2.4 Emergence and Cancellation 

Emergence and cancellation are defined as a pair of flux increase and decrease events of 
different polarities respectively. One can pick up a partial cancellation, where not all but 
only one patch of the involved pair disappears. 

First, we mark a flux change event in time series of flux content in each patch. It is 
defined as a flux increase or decrease which has a duration and a flux change rate more 
than certain threshold. Next, we exclude change of flux content by merging and splitting. 



Cancellations in network field have a time scale of a few ten minutes ( Chae et al. , 2004 



lida et al. , 2010). On the other, flux increase has a shorter time scale, less than 10 



minutes, in emergence event (Otsuji et al. , 2007). We want to set as short threshold for 



the duration as possible to pick up emergences. We set the threshold for duration as 5 
minutes for data set 1 and 10 minutes for data set 2. The threshold for flux change rate 
is taken from a typical flux change rate of cancellation. It is because a cancellation is 
thought to be a more moderate event in flux change rate than an emergence. There are 



some reports of cancellations in network field (Chae et al. , 2002, Park et al. , 2009). Chae 



et al. (2002) reports a flux change rate of cancellation as 3.5 x 10^^ Mx hr ^ from their 



analysis of magnetograms obtained by SOHO /MDl. We use 10^^ Mx hr ^ as a threshold 
for flux change events. The flux change events are marked with these thresholds after 



averaging magnetic flux content over three data point for smoothing. Figure 2.11 shows 
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Figure 2.10: Examples of detected surface processes of magnetic patches, namely, (a) 
merging, (b) splitting, (c) emergence, and (d) cancellation. The background shows the 
magnetic flux density obtained by Hinode/NFl. The red (blue) contours indicate positive 
(negative) patches detected with our threshold. The field of view is 7.2" x 14.4" for all 
the images for all images. 
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Figure 2.11: An example of flux time series contained in one positive patch in data set 1 
and flux change events there. The dashed hne indicates it. The solid hne indicates same 
plot but after removing flux change by mergings and splittings. 



an example of time evolution of flux content in one positive patch in data set 1. Dotted 
line presents time evolution of flux content in the patch. Solid line is the same one but 
after removing flux changes by merging and splitting. Three flux change events, which 
are indicated by arrows, are detected in this patch. 

Then the flux change events in some distance threshold are paired. We set the distance 
threshold for cancellation same as patch tracking, namely 1 pixel in data set 1 and 5 pixels 
in data set 2, because typical proper velocity of canceling patches is not so rapid. On the 
other hand, there is a signiflcant large velocity in emergence. We employ 4 km s~^, which 
corresponds to 4 pixel in data set 1 and 20 pixel in data set 2. 

Our method is capable of detecting partial cancellations, i.e. those in which one or 
both of the involving patches do not disappear. It is also capable for the partial emergence. 
This method categorizes a hybrid event, in which the cancellation and emergence simulta- 
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neously occur, as a single cancellation, a single emergence or a mix of uni-polar processes. 
Thus, the detected number for the emergence and cancellation may be underestimated. 



Figure 2.10 c) shows an example of the emergences. There is a pre-existing positive 
patch in the north of the emerging region. The negative polarity is small and unrecognized 
with our detection condition at 1:05— 1:07UT and grows large enough to be detected at 
1:13UT. The flux content of the positive polarity continues to increase at 1:13UT, which 
is defined as a flux increase in this duration. The distance between patches continues to 
become larger at 1:19— 1:23UT. 

Figure 2.10[ d) shows an example of the cancellations. One negative patch is converging 
to the larger positive patch that is located at the conjunction point of network. They 
contact at 0:49UT and begin to cancel with each other at 0:52UT. The negative patch 
continues to shrink and it becomes smaller than our size threshold at 1:05UT. It totally 
disappears at 1:17UT. On the other hand, positive patch remains after the cancellation 
because the flux content of the positive patch is larger than that of the negative patch. 
So this event is a partial cancellation. 
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Chapter 3 

Statistical Investigations of Magnetic 
Patch Characters 



We investigate statistical properties of magnetic patches in this chapter. Total number 



of patches and total flux amount are investigated in Section 3.1 The flux dependence of 



patch size is investigated to determine the statistical relationship between flux content and 



spatial scale in Section 3.2 Frequency distributions of flux content in both data sets are 



investigated in Section 3.3 We investigate proper velocity of patches, which is necessary 



to evaluate collision frequency of patches, in Section 3.4 Lifetime i.e., a duration for 



which a patch has a flux content larger than our threshold, is investigated in Section 3^ 
We will discuss these results and show our interpretations in Section 5.2 and 5.3. 



3.1 Total Number and Total Flux Amount of Patches 

Total number and total flux amount of patches are investigated. Table 4.1 shows the 
results of both data sets. 

26745 positive and 26459 negative are detected in data set 1. They correspond to 
~ 130 patches of one polarity in one image. As for data set 2, 97018 positive and 92102 
negative patches are found in data set 2. These numbers are sufficient for statistical 
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Data Set 1 


Data Set 2 










26745 


97018 


Positive 


(134.4 per frame) 


(59.1 per frame) 




26459 
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(133.0 per frame) 


(56.1 per frame) 
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Tracked Patches 






Positive 


loob 


zlozo 


Negative 


1637 


19544 


Total 


3273 


41367 


Total Flux Amount 








i.DO X iU iViX 


i. ( U X iU iVlX 


Positive 


(2.53 Mx cm-2) 


(3.76 Mx cm-2) 




2.32 X 10^0 Mx 


1.54 X 10^0 Mx 


Negative 


(3.60 Mx cm-2) 


(3.42 Mx cm-2) 




3.95 X 10^0 Mx 


3.42 X 10^0 Mx 


Total 


(6.13 Mx cm-2) 


(7.16 Mx cm-2) 



Table 3.1: Number of patches and flux amount in the data sets. 



analysis. They correspond to ~ 60 patches of one polarity in one image. 1636 positive 
and 1637 negative patches are detected as one tracked patch in data set 1. 21823 positive 
and 19544 negative patches are detected as one tracked patch in data set 2. The numbers 
of tracked patches may also be sufficient for statistical analysis. 



Figure |3.1| shows time series of total flux amount in the observational period in the 
target regions. The red lines denote total flux amount of positive polarity. The blue ones 
denote those of negative polarity. The black ones denote sum of them. Time series of 
total flux amount in data set 1 is shown in the upper panel. The averaged total flux in 
data set 1 is 1.63 x 10^° Mx, 2.32 x 10^° Mx and 3.95 x 10^° Mx for positive, negative, and 
both polarities, respectively. They correspond to 2.53 Mx cm~^, 3.60 Mx cm~^ and 6.13 
Mx cm~^ as flux density averaged over whole area. Negative flux content is larger than 
that of positive one by 42% in this region. The maximum differences from average are 
8.3% and 5.1% for positive and negative polarities. Total flux amounts of both polarities 
are constant with an accuracy of 10%. Same plot for data set 2 is shown in the lower 
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panel. There are many clear large flux emergences (but not so large as that forming active 
region) in data set 2. The averaged total flux content is 1.70 x 10^° Mx, 1.54 x 10^° Mx 
and 3.42 x 10^° Mx for positive, negative, and both polarities. They correspond to 3.76 
Mx cm~^, 3.42 Mx cm~^ and 7.16 Mx cm~^. The difference between polarities is 9% The 
maximum differences from average are 45% and 52% for positive and negative polarities. 



3.2 Dependence of Patch Area on Flux Content 



Figure 3.2 shows a scatter plot of flux content and area of patches in data set 1. There is 



a positive correlation between them. We fit with a power-law function, namely 

s(0)=so f^y. (3.1) 



where 5*0 is the reference area, 0o(= 10^^ Mx) is the reference magnetic flux content, 
and a? is the power-law index. We obtained 5*0 = (3.33 ± 0.12) x 10^® cm^ and a? = 



0.711 ±0.001 with the la-error of fitting. Solid line in Figure 3.2 denotes the fitting result 



3.3 Frequency Distribution of Flux Content 



Figure 3^ shows a frequency distribution of the flux content in data set 1. The red and 
blue dashed histograms denote frequency distributions of positive and negative polarities. 
The solid histogram denotes the sum of them. The distribution drops down below </)th = 



10^''-^ Mx, which is suggested as a detection limit defined in Section 5.1.3 Only a limited 
number of the patches are detected below it. The vertical green line denotes it. We set 
0th = 10^'^'^ Mx as a limit of statistical analysis in data set 1 from here. We make a fitting 
of the power-law function between 10^^'^ Mx and 10^^ Mx with a least-squares fitting, 
namely 

n(0) = no (-^) ' (3.2) 
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Figure 3.1: Time series of the total flux amount contained in the field of view during the 
observational period in data set 1 (upper) and data set 2 (lower). The red and blue line 
denote total flux amount of positive and negative polarities, respectively. The black line 
denotes sum of them. 
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3.3. FREQUENCY DISTRIBUTION OF FLUX CONTENT 



Flux Content vs Size 




i20 



Flux Content [Mx] 



Figure 3.2: Scatter plot of the flux content and the area of magnetic patches in data set 
1. The red (blue) asterisks denote those of positive (negative) polarity. 200 data points 
for each polarity are randomly selected for the plot. The solid line denotes the fltting 
result with an assumption of a power-law distribution. 
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where Uq is the reference frequency density, 0o(= 10^^ Mx) the typical scale of reference 
magnetic flux content, and 7 is negative value of a power-law index of the distribution. 
The black dashed line indicates the fitting result. We obtain no = 2.42 x 10"'^^ Mx~^ 



cm~ and 7 = 1.78 ± 0.05 with the la-error of fitting. Parnell et al. (2009) reports that 
error of a power-law index from the selection of fitting range is larger than that of fitting 
itself. We evaluate that error by changing the minimum value of fitting range. The index 
varies between —1.78 and —1.91 with the minimum value from 10^'''^ Mx to 10^* Mx. the 
total error of our fitting of the power-law index becomes 0.18 with taking it into account. 



The same analysis is done for data set 2. Figure 3.4 shows a frequency distribution 
in data set 2. The distribution is dropping down below 0th = 10^"^'^ Mx. The difference 
of the detection limit may be caused by the fact that there are some emerging fluxes in 
this region during the observational period. We obtain Uq = 1.63 x 10^'^^ Mx^^ cm~^ and 
7 = 1.93 ± 0.03 with the la-error of fitting. The index varies between —1.90 and— 1.94 
with a minimum fitting range of 10^^'^ Mx to 10^® Mx. The total error of our fitting of 
the power-law index becomes 0.07. 

These results are consistent with Parnell et al. (2009), who report uq = 1.5 x 10"'^^ 
Mx~^ cm"^ and 7 = —1.85 ± 0.14 with 0o = 10^^ Mx for one polarity. 



3.4 Proper Velocity 

We investigate proper velocity of magnetic patches. The aim here is to investigate average 
velocity of patches in our data set. The proper velocity is important to evaluate the 



collision frequency of patches (Schrijver et al. , 1998). Proper velocity is defined as a 



velocity of the center of a patch in this study. The center of patches is determined by 
weighting pixel position with magnetic flux, namely 



(3.3) 
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Figure 3.3: Frequency distribution of the magnetic flux content in data set 1. The red 
(blue) dashed histogram denotes that of positive (negative) polarity. The solid histogram 
denotes the sum of them. The black dashed line indicates the fitting result of it. The 
obtained power-law index is —1.78 with a fitting range between lO^'''^ Mx and 10^^ Mx. 
The vertical green line denotes the detection limit, 0th = 10^^'^ Mx. 
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Frequency Distribution of Flux Content 
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Figure 3.4: Same plot as Figure 3^ but for data set 2. The power-law index is —1.90 
with a fitting range between 10^^ '' Mx and 10^^ Mx. The vertical green line denotes the 
detection limit, (hfh = 10^^"^ Mx. 
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Here xl is a position of center of i-th patch at j-th time step, N- is a pixel number of i-th 
patch at j-th time step, B^{x) is a flux density at the position of x at j-th time step, 
and cc^^ is a position of A;-th pixel of i-th patch at j-th time step. The proper velocity of 
patches at merging or splitting removed in the following statistical analysis since, in such 
timing, an artificially large displacement may occur. 



Figure 3^ shows a number distribution of proper velocity. The red dashed line in- 
dicates that of positive patches. The blue one indicates that of negative patches. The 
black solid one devotes sum of them. They reach maximum value around 10^ cm s~^. 
There is an enhancement in the range of high velocity, namely from 10^ cm s~^ and 10® 
cm s~^. We found 123 and 124 points out of 24096 positive and 24068 negative points 
there. They correspond to ~0.5% of data points. Surface processes of patches involving 
undetected patches, which we cannot remove in this analysis, may cause them. Then we 
remove proper velocity larger than 10^ cm s~^ and evaluate average proper velocity of 
patches. There are also enhancements between 3 x 10^ and 2 x 10^ cm s^^. Such rela- 
tively rapid motion of patches might corresponds to a peculiar events such as emergence 
driven by buoyancy, for which further investigation is necessary. The obtained averaged 
proper velocities of positive, negative, and both patches are 1.10 km s~^, 1.19 km s~^, and 
1.14 km s~^, respectively. The medians of proper velocity are 0.81 km s~^, 0.79 km s~^, 
and 0.80 km s~^ for positive, negative, and both patches, respectively. They are slightly 
smaller than the averages due to the peak offset in the distribution to the larger velocity 
side. 



Further a relationship between proper velocity and flux content is investigated. Figure 



3.6 shows a scatter plot between them. The red and blue asterisks denote data point of 



positive and negative patches. The green horizontal line denotes the threshold to cut the 
large velocity enhancement. 300 and 8 data points for each polarity are randomly selected 
for plot below 10'' cm s^^ and above it, respectively. We make a least-squares fitting with 
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a power-law function. The fitting form is written as 



^(</')=^o(^^J (3.4) 

where Vq is the reference velocity, (po is the reference magnetic fiux which we set 10^^ Mx, 
and is the power-law index of the dependence. We obtain Vq = (7.14 ± 1.12) x 10^ 
cm s~^ and = —0.227 ± 0.003. The fitted result is shown by the black dashed line in 
Figure [3^ 



3.5 Frequency Distribution of Lifetime 

We investigate the frequency distribution of lifetime of patches in this section. We define 
lifetime as a duration for which a patch is recognized based on our detection, namely from 
the appearance to the disappearance in two-leveled images. In data set 1, the averaged 
lifetime of positive, negative and both polarities are 17.4 minutes, 17.2 minutes, and 17.3 



minutes, respectively. Figure 3.7 denotes a number distribution of lifetime of tracked 
patches. The red, blue, and black lines denote those of positive, negative, and both 
polarities. We make a least-squares fitting of a power-law function in the range between 
6 minutes and 30 minutes. The fitting form is 



T 



N{r) = No . (3.5) 

where Nq is the reference number, Tq is the reference lifetime that we set 10 minutes, and 
is the power-law index of the dependence. We obtained Nq = (6.51 ± 0.05) x 10 and 



= —1.45 ± 0.13 with a Icr-error of fitting. The dashed line in Figure 3.7 represents 
the fitting result for the distribution of both polarities. The power-law index steeper than 
— 1 in number distribution means that the lifetime of patches becomes shorter as their 
smaller fiux content. 
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Figure 3.5: Number distribution of the proper velocity. The red (blue) dashed line 
indicates that of positive (negative) patches. The black solid line indicates sum of them. 
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Figure 3.6: Scatter plot of the flux content of the patches and their proper velocity in 
data set 1. The red (blue) asterisks denote data point of positive (negative) patches. The 
green horizontal line denotes the threshold to cut the large velocity enhancement. The 
black dashed line indicates a fitting result with a power-law function. 
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3.5. FREQUENCY DISTRIBUTION OF LIFETIME 



Further dependence of lifetime on flux content is investigated. We define flux content 



of tracked patches by averaging flux content over whole their lifetime. Figure |3.8| shows 
the result of it. Red and blue dashed lines indicate those of positive and negative polar- 
ities, respectively. Solid line indicates that of both polarities and horizontal dashed line 
indicates the observational period of this data set. Vertical dashed lines denote detection 
limit in this study, namely (ptt = 10^^'^ Mx, and twice of it. 



There are two domains from the discussion in Section [5. 1.3[ One is below 2<^th, where 
we cannot distinct the apparent disappearance from the actual disappearance, and the 
other is above it, where one apparent disappearance cannot remove the patch. There 
seems to be a change of dependence around it. However, the lifetime hits near the ob- 
servational duration around 10^^ Mx. To investigate this effect, we investigate a scatter 



plot of magnetic flux content of tracked patches and their lifetime. Figure |3.9| shows the 
result. Red and blue asterisks indicate data points of positive and negative polarity. Total 
numbers of positive and negative patches are 1636 and 1637 respectively. The horizontal 
dashed line presents the whole observational period of this data set. Some of data points 
hit the observational duration around 10^^ Mx. It is difficult to distinct the effect of it 
and physical change of dependence. 

Same analysis is done for data set 2. The observational period is much larger than 
that of data set 1. The averaged lifetime is 22.5 minutes, 23.9 minutes, and 23.2 minutes 



for positive, negative, and both polarities. Figure 3.10 shows a number distribution of 
lifetime. It seems a power-law distribution in the range below 70 minutes. We make a 
least-squares fitting in the range between 30 — 60 minutes. Nq = (1.34 ± 0.21) x 10^ and 
= —1.53 ± 0.06 are obtained with tq = 40 minutes. There is a rapid drop above it. 



which is different from data set 1. Figure 3.12 shows a scatter plot of flux content and 



lifetime. The horizontal dashed line shows the observational duration of this data set. We 



see that the longest lifetime is shorter than the observational period. Figure 3.11 shows a 



dependence of averaged lifetime on flux content. There is a change of slope around 20th- 
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Figure 3.7: Number distribution of the lifetime of the detected patches in data set 1 
in log-log plot. The red/blue/black line denotes that of positive/negative/both polarity. 
The dashed line indicates a fitting result of both polarities. It has a power-law index of 
-1.45. 
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Figure 3.8: Dependence of the averaged lifetime on the magnetic flux content in data set 
1. The red and blue lines indicate those of positive and negative polarities respectively. 
The solid line indicates that of both polarities. The horizontal dashed line indicates the 
observational period of this data set. The vertical dashed lines denote 0th and 20th in 
data set 1. 
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Figure 3.9: Scatter plot of the magnetic flux content and the lifetime in data set 1. The 
red (blue) asterisks indicate data points of positive (negative) polarities. The horizontal 
dashed line indicates the whole observational duration of this data set. 
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3.5. FREQUENCY DISTRIBUTION OF LIFETIME 



The lifetime decreases again below it. However, the lifetime becomes nearly independent 
of flux content above it. The value of constant lifetime is ~ 60 minutes, which is similar 



value of upper limit of power-law distribution in Figure 3.10 
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Figure 3.10: Same plot as Figure 3^ but for data set 2. The power-law index of the 
fitting line is —1.53. 
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Figure 3.11: Same plot as Figure 3.8 but for data set 2 
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Figure 3.12: Same plot as Figure 



3.9 



but for data set 2. 
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Chapter 4 



Frequencies and Flux Amounts of 
Surface Processes of Magnetic Patches 



Frequencies and flux amounts of surface processes are investigated in this chapter. Total 
frequencies and flux change rates of processes in both data sets are summarized in Section 



4.1 The purpose of this investigation is to determine what processes is dominant in 



surface processes. Further, flux dependences of merging, splitting, and cancellation are 



investigated in Section 4.2 We will construct a physical picture of flux maintenance in 



quiet regions from these results in Section 5.3 



4.1 Total Frequencies and Flux Amounts of Surface 
Processes 

We investigate changing rates of total flux amount by four surface processes of magnetic 
patches. Along with the total numbers of patches, their flux density averaged over the 
field of view and the observational period (in Mx cm~^), the frequencies of events in the 



magnetic flux density (in Mx cm ^ s ^) are summarized in Table 4.1 



The total numbers of merging for positive and negative patches are 493 and 482 in 
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data set 1. They are 5905 and 4252 in data set 2. The total number of sphtting has a 
same order of magnitude. Those of sphtting for positive and negative patches are 536 and 
535 in data set 1. Those in data set 2 are 5764 and 4143. The total numbers of emergence 
in data set 1 and in data set 2 are 3 and 27 respectively. They are much smaller than 
those of merging and splitting in both data sets. The total numbers of cancellation are 
86 and 775, which are again smaller than those of merging and splitting. 



To determine what processes dominate a flux balance in the region, we investigate 
change rates of total flux amount by each process. The occurrence time is obtained by 
dividing the total flux by this rate. Note that these are averages over the flux content 
in patches. The flux dependence of the occurrence times is discussed later in Section 



4.2 The results are shown in the parenthesis in Table 4.1. The change rates of flux 



density involving merging process are deflned as sum of flux content of patches before 
the events. Those in data set 1 are 1.65 x 10^'^ Mx cm^^ s^^ and 3.53 x 10^"^ Mx 
cm~^ s^^ for positive and negative patches, respectively. Time scales of those processes, 
^tot/id^tot/dt\^^^), are evaluated as 1.53 x 10^ s and 1.02 x 10^ s from these values. We 
obtain 3.36 x 10^'^ Mx cm^^ s^^ and 1.77 x 10^^ Mx cm~^ s^^ for positive and negative 
patches in data set 2. They correspond to the time scales of 1.21 x 10'^ s and 1.93 x 10^ 
s, respectively. The time scales of these processes are much shorter than the estimated 



time scales by cancellation and emergence reported in previous studies (Livi et al. , 1985 



Schrijver et al. 1998 Hagenaar 2001 Thornton & Parnell, 2011). The change rates of 



flux density involving splitting process are deflned as sum of flux content patches before 
the events. We obtain 1.48 x 10~^ Mx cm"^ s~^ and 3.03 x 10"^ Mx cm"^ s"^ for positive 
and negative patches in data set 1, which corresponds 1.71 x 10'^ s and 1.19 x 10'^ s, 
respectively. In data set 2, we obtain 2.38 x 10^^ Mx cm"^ s"^ and 1.72 x 10"^ Mx cm~^ 
s~^ for positive and negative patches. The time scales of positive and negative patches 
are 1.58 x 10^ s and 1.99 x 10^ s. They are same order of those of merging in both data set 
and much shorter than those of emergence and cancellation in the previous studies again. 
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57 PROCESSES 





Data Set 1 
Positive Negative 


Data Set 2 
Positive Negative 


Tracked Patch 


1636 


1637 


21823 


19544 




(2.53) 


(3.60) 


(3.76) 


(3.42) 


Merging 


OoO 


OoO 


5905 


4252 




(1.65 X 10-3) 


(3.53 X 10-3) 


(3.36 X 10-3) 


(1.77 X 10-3) 


SnlittinB" 


493 


482 


5764 


4143 




(1.48 X 10-3) 


(3.03 X 10-3) 


(2.38 X 10-3) 


(1.72 X 10-3) 


Emergence 




Q 

o 




27 


C a n CP ] ] a t i on 




86 


775 








(2.09 


X 10-5) 


Unipolar Events 










Increase& 










Appearance 


503 


556 


6185 


7633 


Decrease& 










Disappearance 


381 


398 


5440 


5400 



Table 4.1: Summary of total numbers of the patches and surface processes in our data 
sets. Flux density and change rates of the flux amount are shown in the parentheses with 
units of Mx cm"^ and Mx cm"^ s-^, respectively. 



The timescale of cancellation is investigated in data set 2. We detect 610 cancellations. 
Change rate of flux amount is 2.09 x 10"^ Mx cm-^ s-^, which corresponds 1.71 x 10^ s 
as a time scale of cancellation. Time scales of emergence and cancellation are in the same 
order with the previous studies. It is much longer than those of merging and splitting. 
These results deduce that merging and sphtting are dominant surface processes in quiet 
regions. 

We should note that there are also many unipolar increases and decreases in both 
data set. Total numbers of the detected unipolar processes are shown in the last two 
hues in Table 5.1. These events are interpreted as the surface processes with the patches 
below the detection limit. We show some examples of these processes and discussions in 
Appendix B to support the interpretation. 
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4.2 Frequency Distributions of Interactions of Magnetic 
Patches 

4.2.1 Merging 

We investigate flux dependences of occurrence of mergings in data set 1 namely, the 
frequency distribution and the probability distribution. The reason using data set 1 is 
that the time interval is short enough compared with the occurrence time (~17 minutes). 

First, we investigate the frequency distribution of mergings. The freqeuncy distribu- 
tion of merging is deduced by 

n„APP yytot 

mrg _ mrg (A 

dt ~ t'-'SAf ^ ■ ' 

where iV^rg denotes the total occurrence of mergings throughout the data, t*°* the obser- 
vational period, S the observational area, and A0 the bin size in the flux content. The 



result is shown in Figure 4. 1 The distribution seems to have a power-law dependence on 



flux content and can be a least-squares fitted in the forms as, 

'^o.mrg is the reference frequency at (po, 0o(= 10^^ Mx) is the reference flux content, and 
Q^mrg is the power-law index of the probability distribution of merging. We obtained 
^^o,mrg = (8.14 ± 0.86) X 10-3^ Mx~i cm-2 and a^rg = 1-66 ± 0.07 for the positive 
patches, no,mrg = (7.41. ± 0.72) x 10"^^ Mx^^ cm^^ g-i ^^^^ ^ i ^O ± 0.06 for the 
negative patches, and no.mrg = (1.60±0.77) x 10^'^'' Mx"-*^ cm^^ s^^ and ttmrg = 1.68±0.05 
for both patches. The errors mean a Icr-error of the least-squares fitting. 

Next, we investigate the probability distribution of merging. The apparent probability 
distribution of merging is obtained by 

mrg _ -L ^ "mrg / . „n 

dt ~ n(0) dt ■ ^ ^ 
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Flux Content Before Merging [Mx] 

Figure 4.1: Apparent frequency distribution of mergings with 0th = 10^^'^ Mx in data set 
1. It is made from merging where the parent patches have a flux content of more than 0th- 
The red/blue/black solid lines indicate observational results for positive/negative/both 
patches. The red/blue/black dashed lines indicate fitting results with a range of lO^^'^- 
— 10^^ Mx. The power-law indexes of the fitting lines are —1.66., —1.70, and —1.68. 
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Figure |4.2| shows the resuh. We make a least-squares fitting in a range of 10^^'^ — 10^^ 



Mx, where the number of the merging event is enough for the fitting. The fitting form is 



dP, 



APP / , \ /3n„. 



mre 



^0,mrg ^ (4.4) 



where Po,mrg is the reference probabihty at (pQ, 0o(= 10^^ Mx) the reference flux content, 



and /3inrg the power-law index. We obtained Po,mrg = (2.52 ± 0.08) x 10 s and /3 



mrg 



0.28±0.05 for the positive patches, Po,mrg = (2.52±0.08) x 10"^ s^^ and /3mrg = 0.26±0.05 
for the negative patches, and Po,mrg = (5.12 ± 0.11) x 10^"^ s^^ and /3mrg = 0.28 ± 0.04 for 
both patches with (pQ = 10^® Mx. The errors mean la-errors of least-squares fitting. 



4.2.2 Splitting 

We investigate the flux dependence of splitting occurrence in the same manner as that 
for mergings in the previous section. 

First, we investigate the frequency distribution of splittings. It is 

_ A^pt .... 
dt P^^SAf ^ ■ ^ 



The result is shown in Figure 4.3 The distribution seems to be a power-law dependence 



on flux content again and we make a least-squares fitting in the forms as, 

-- "-o.spit — • (4.6) 



dt 



'Jo 



We obtained no,mrg = (3.24 ±0.62) x 10"^^ Mx^^ cm^^ g-^ and a^rg = 1.01 ±0.07 for the 
positive patches, no,mrg = (3.41 ± 0.55) x 10"^^ Mx"-"^ cm"^ s"-*^ and a^^-g = 1.14 ± 0.08 
for the negative patches, and r;,o,mrg = (6.75 ± 0.95) x 10~^^ Mx~^ cm~^ s^^ and amrg = 
1.01 ± 0.07 for both patches. 

Next, we investigate an apparent probability distribution of splitting in data set 1. 
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Figure 4.2: Apparent probability distribution of mergings with 0th = 10^'''^ Mx in data set 
1. It is made from merging where the parent patches have a flux content of more than 0th- 
The red/blue/black solid lines indicate observational results for positive/negative/both 
patches. The red/blue/black dashed lines indicate fitting results with a range of 10^'''^ 
— 10^^ Mx. The power-law indexes of the fitting lines are 0.28, 0.26, and 0.28. The 
horizontal and vertical dashed lines indicate time scale of 33 minutes and the detection 
limit, respectively. 
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Flux Content Before Splitting [Mx] 

Figure 4.3: Apparent frequency distribution of splittings with 0th = lO^"^'^ Mx in data set 
1. It is made from splitting where the parent patches have a flux content of more than 0th- 
The red/blue/black solid lines indicate observational results for positive/negative/both 
patches. The red/blue/black dashed lines indicate fitting results with a range of lO^'''^- 
— 10^^ Mx. The power-law indexes of the fitting lines are —0.89, —1.14, and —1.01. 
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The apparent probability distribution of splitting is obtained in the same manner as that 
of merging. The apparent probability distribution of splitting is induced as 

dt n{cf)) dt ' ^ ■ ^ 



Figure 4A shows the result. The strong increase in the range larger than 10^^ Mx is 
probably caused by the lack of a patch number in the analysis. On the other hand, there 
is a drop in the flux range near 0th where the number of patches is enough for a statistical 
study. We interpret this dropping as an effect of splitting into the area below ^th- We 



quantify this effect in the discussion in Section |5.3.2[ We see that the probability of 
splitting is almost constant as 1.0 x 10"'^ s~^, which means a time scale of 33 minutes, 
in the range enough above 0th) 3.0 x 10^^ — 1.0 x 10^^ Mx. It means that frequency of 
splitting is independent of the parents' flux content. 



4.2.3 Cancellation 

We investigate an apparent frequency distribution of cancellations in this section. Because 
the total number of cancellations in data set 1 is too small, we use data set 2. The 
frequency distribution of cancellation is obtained in the same manner as the merging and 
splitting, namely 

n„APP Tytot 

cnc cnc q\ 

dt ~ tt°tSA0' ^ ■ ' 



Figure 4.5 shows the result. We make a least-squares fitting with a form of 

^^^nc^ _ dnp^cm 

dt ~ dt 



(4.9) 



We obtained dno^^^Jdt = (1.29 ±0.12) x 10"^^ Mx'^ cm'^ s'^ and 7cnc = 2.48 ±0.26 and 
the fitting range from 10^^''' Mx to 10^*'^ Mx. The errors indicate Icr-errors of least-squares 
fitting. 
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Figure 4.4: Apparent probability distribution of splittings with (/)th = lO^'''^ Mx 
in data set 1. The red/blue/black solid lines indicate observational results for posi- 
tive/negative/both patches. The blue and black dashed lines indicate analytical curves 
with ko = 5.0 X 10~^ s~^ and 1.0 x 10~^ s~^, which are discussed in Section 5.3.2 The 
horizontal dashed line indicates a time scale of 33 minutes. 
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Figure 4.5: Frequency distribution of cancellations on decreased flux in data 2. The 
histogram and dashed line indicate the observational result and the fitting result with 
a power-law distribution, respectively. The vertical green line shows detection limit, 
0th = lO"'^'^'^ Mx. The fitted power-law index is —2.48. 



CHAPTER 4. FREQEUNCIES OF SURFACE PROCESSES 66 



We investigate the dependence of the power-law index on the bin size for the least- 
squares fitting because the total number of cancellation is not so large in data set 2. The 



result is shown in Figure 4^ The error bars indicate Icr-errors of the least-squares fitting 
in each case. We see that the error has the tendency to increase with the larger bin 
size, which is caused by the decrease of the fitting points. The proper power-law index is 
roughly from —2.3 to —2.6. 

Then, we investigate an apparent probability distribution of cancellation. It is given 

as 



cnc _ ' ^cnc /a i q\ 

dt ~ n(0) dt ' ^ ^ ^ 

Figure 4.7 shows the result. We make a least-squares fitting in a range of 10^^'^ — 10^^'^ 



Mx in the same manner as that for merging and splitting. The fitting form is 



npAPP 

-P0,cnc(^) (4.11) 



dt 



'Jo 



We obtained Po.cnc = (1-52 ± 0.15) x 10"^ s^^ and /5cnc = -0.50 ± 0.24. The error bars 
mean la-errors of the least-squares fitting. 
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Figure 4.6: Dependence on fitting bin size of the power-law index of cancellation in data 
2. The error bars show Icr-error of the fitting. 
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Figure 4.7: Probability distribution of cancellation in data set 2. The histogram and 
dashed line indicate the observational result and the fitting result with a power-law distri- 
bution, respectively. The vertical green line shows detection limit, 0th = 10^^'^ Mx. The 
fitted power-law index is —0.50. 
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Chapter 5 



Discussions and Interpretations 



We discuss the results of Chapter 3 and Chapter 4 in this chapter. First we summarize 
and expand magneto-chemistry equation (M-C equation hereafter), which is suggested 



by Schrijver et al. (1998) and useful to evaluate the observational results, in Section 5.1 



Some forms connecting observational results and elementary functions in M-C equation 
are deduced. Discussions of characteristics of patches and frequencies of surface processes 



are done in Section [5^ and [573| We suggest a physical picture of flux maintenance on the 
solar surface from these discussions in Section 15.41 



5.1 Magneto-Chemistry Equation 

The detailed description of magneto-chemistry equation is summarized and expanded 
in this chapter. We deduce some forms connecting frequency distributions of surface 
processes with observables, namely total number of patches, total flux amount, total 
numbers of surface processes, time scales by surface processes, and averaged lifetime. A 



general description of magneto-chemistry equation is summarized in section 5.1.1 We 



deduce general forms in section [5.1.2 In section 5.1.3, we take the detection limit of 



analysis into the account. We stand on one interpretation of unipolar processes that they 
are surface processes between magnetic patches below the detection limit. 
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5.1.1 General Description 



The magneto-chemistry equation was proposed by Schrijver et al. (1997), which describes 
a relation among frequency distribution of flux content, n{(j)), and those of four surface 
processes. This equation treats magnetic structures as isolated patches. This treatment 
is valid with continuity and no dissipation of the field, here B. The remained condition, 
no dissipation of magnetic field, is probably valid because of the large magnetic Raynolds 
number, Rm = LoVo/r], on the soar surface. We obtain Rm ~ 10^ — 10'' in the solar 



photosphere, assuming partial ionization and collision dominated plasma (Spitzer, 1962 



Priest, 1984). Then magneto-chemistry equation is written as 

dn±{(j)) 



dt 



where n±{(f)) denotes a frequency distribution of positive (negative) patches with a flux 
content (f), S±{(f)), L±{(j)), K±{(j)), and M±{(f)) denotes changes of number density of pos- 
itive (negative) patches with a flux content of by emergence, merging, splitting, and 
cancellation, respectively. We describe them with probability density distributions of ac- 
tivities, l{x,y), k{x,y), and m(x,?/). Here, l{x,y) is defined as a probability density of 
mergings where patches with flux contents of x and y merge into a patch with a flux 
content of x + y, k{x, y) is defined as a probability density of splittings where a patch 
with flux content x + y splits into patches with flux contents of x and y, and m{x,y) is 
defined as a probability density of cancellations where patches with flux contents of x and 
y cancel and results in a patch with a flux content of |a;-|-?/|. By this definition, one obtain 

1 /"^ 

L± = - n±{x)n±{(j) — x)l{x,(j) — x)dx 



2 







n±{(j)) I n±{x)l{(f), x)dx (5.2) 




/•oo 

K± = 2 n±{x)k{(j),x — (j))dx — n±{(j)) I k{x,(f) — x)dx (5.3) 
Jo Jo 
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M4 



n± (0 + x)nzp {x)m{(j) + x, x)dx 



—n±{(f)) / n^{x)m[(j),x)dx. 
'0 



(5.4) 



The first and second terms in equations (5.2)-(5.4) are production and loss terms respec- 



tively. We remove difference between polarities as 



n+(0 


) = ^-(0) 


= n{(p) 




) = 


= S{<P) 




) = 1-{<P) = 


-m 




) = k4<f>)- 


= m 


m+(0 


) = m_(0) 


= m{(f)) 



(5.5) 
(5.6) 
(5.7) 
(5.8) 
(5.9) 



Equation (5.5) is applicable only in well mixed-polarity and flux balanced region, which 



is probably apphcable in quiet regions. 

5.1.2 Relationships with Observables 

We show specific forms connecting probability density distributions, namely l{x, y), k{x, y), 
and m{x,y), and observables in this section. These forms are used in the discussion of 
chapter 4 and 5. 

First, we summarize the forms of total number of patches, A^tot, and total flux amount. 



$tot; in the system from Schrijver et al. ( 1997). We assume the homogeneity of a frequency 



distribution in this study. Then we obtain 



A^tot = / n{(j))d(f) 
Jo 

POD 

$tot = / (l)n{(f))d(p. 
Jo 



(5.10) 
(5.11) 



These variables are more directly comparable with the observables. 
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Next we deduce forms representing probability density distributions, total occurrence 
rates, and total numbers of surface processes. Probability density distribution means 
probability density of one patch. It would represent the physical character more clearly. 
Frequency distributions of merging and splitting in one polarity equal to loss terms in 



Equation (5.2) and Equation (5.3). On the other hand, that for cancellation equals twice 
of loss term because we have to sum up occurrence of cancellation in both polarity. We 
obtain them as 



^"""^■■§.(0) = n{(f)) I n{x)l{(j),x)dx (5.12 



dt 

(0) = n(0) / k{x,(j) — x)dx (5.13) 







dt 

dricnc 

™^ (0) = 2n(0) / n{x)m{(f),x)dx (5-14) 
dt Jo 

where dn^^j-g/dt, duspu/dt and dricnc/ dt are change rates of frequency distribution by 
merging, splitting, and cancellation respectively. Frequency distribution of emergence, 
duemrg/dt, cquals to a source term itself from their definitions, namely 

— ^^(0) = 2n(0) / n{(f) + x)m{(f),(f) + x)dx. (5.15) 
dt Jo 

We obtain probability density distributions simply by dividing frequency distributions by 
frequency distribution of flux content as 



^^(0) = -l-^!^ = 2 rn{x)n{<f) + x)m{<P,<l) + x)dx (5.16) 
dt ~ 

^(0) = I nix)li<P,x)dx (5.17) 



1 




n(0) 


dt 


1 


dfl^arg 


n(0) 


dt 


1 


drispit 


n{(f)) 


dt 


1 


dflcnc 


n{(f)) 


dt 



oo 



(0) = rTTV^sF = / Hx,(f)-x)dx (5.18) 
Q^'^ (0) — ZTTV^tT"" ~ ^ / n{x)'m{(l),x)dx (5.19) 

where dPemrg/dt, dPrarg/dt, dPspit/dt and dPcnc/dt are probability density distributions of 
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emergence, merging, splitting and cancellation respectively. We evaluate total occurrence 



rates of them. By integrating Equation (5.12)-(5.15) on from to oo, we obtain 



dt 



cmrg 



dt 



dt 



mrg 



spit 



n{(f)) I n{(j) + x)m{(j), (j) + x)dx 





A 



n{(j)) / n{x)l{(f)^ x)dx 
Jo 

r<t> 

n{(l)) / k(x,(f) — x)dx 



dt 



n{(f)) I n{x)m{(f),x)dx 





(5.20) 
(5.21) 
(5.22) 
(5.23) 



Further we want to deduce total change rates of flux amount by surface processes. By 



multiplying Equation (5.12)-(5.15) by and integrating on (p from to oo, we obtain 



dt 

9$tot 



PCX) 






0n(0) / 


cmrg 


"'0 



dt 

5$tot 



dt 

a$tot 



POO 








mrg 




POO 








spit J 


"'0 



dt 



4>n{(f)) I n{x)m{(f),x)dx 





(5.24) 
(5.25) 
(5.26) 
(5.27) 



The time scale for each process are defined by dividing total flux amount. Equation 



(5.11), by total change rate of flux amount , Equation (5.24)-(5.27). 



emrg 



tot 



9$ 



tot 



dt 



-\ -1 



(f)n{(j))d(j) 







(pn^cp) I n{(f) + x)m{(l), (p + x)dx 





(5.2J 
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9$ 



n -1 



tot 



tot 



dt 



mrg 



(j)n{(f))d(j) 







Tsplt 



tot 



(j)n{(l)) / n{x)l{(f),x)dx 
Jo 

-\ -1 

^-^-tot 



(5.29) 



9t 



spit 
oo 



(f)n{(j))d(j) 



(j)n{(l)) I k{x,(f) — x)dx 





(5.30) 



P$tot 




9t 


cnc- 





poo 

(pn^cfydcf) 




'J 


POO 




0n(0) j 


n(x)m((f), x)dx 




(i0 



(5.31) 



At the last step, we deduce lifetime of patches. Patches disappear at merging or 
canceling with a patch that has larger flux content. The probability of disappearance is 



written from Equation (5.17) and Equation (5.19) as 



disapp 



dt 



1 f)P f)P 

/ i\ -L L/-I tcnc , '^^ tmre 



n{x)m{(f), x)dx + / n{x)l{(j),x)dx. 



(5.32) 



The lifetime is an inverse of probability of disappearance, namely 



r(0) 



disapp 



dt 



n{x)m{(j)^x)dx + / n{x)l{(f), x)dx 



(5.33) 
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5.1.3 Consideration of Detection Limit 

We want to introduce one idea, detection limit here. As is explained in Chapter 2, there 
are many unipolar appearances and disappearances in the observational data. One possi- 
ble interpretation is surface process among involved patches below the observational limit 



as reported by Lamb et al. (2010). In a statistical analysis, it probably leads a misun- 
derstanding when we investigate in the range where most patches are missed. However, 
we cannot avoid missing the patches with smaller size than spatial resolution and weaker 
signal than observational limit. We introduce an analytical limit of flux content, 0th, to 
evaluate the effect quantitatively. The term a detection limit is defined as a flux content 
above which we can pick up all patches. We obtain it from the frequency distribution of 
flux content. We obtain 0th ~ 10^^ for MDI FD data (Figure 1.3). We cannot justify the 



assumption completely. However, one justification of the interpretation is that frequency 
distribution by SOT/NFI, namely high-resolution data, belongs to the same power-law 
distribution although it is flattened in MDI FD line. 

With a limit of statistical analysis, the forms for observables should be changed. We 
can pick up processes where all involved patches have flux content larger than 0th- The 
number density of patches, flux density, total number of patches and total flux amount 
are obtained simple by changing the range of integration from [0,oo] to [0th, cxd], 

n(0)# (5.34) 

rth 

POO 

= A / 0n(0)rf0. (5.35) 
•J 4>th 

where superscript APP means apparent observable taking into the detection limit into 
account. 

On the other hand, apparent frequency distributions of surface processes become 

^^APP ^oo 

— ^^(0) = 2n(0) / ?2(0 + x)m(0, + x)(ix (5.36) 

Ot In 
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Q APP 

Q APP 



n(0) / n{x)l{(j), x)dx 
r4>-<f>th 

n{(t)) I k{x,(p — x)dx 

poo 

2n{(f)) / n{x)m{(f),x)dx. 



(5.37) 
(5.38) 
(5.39) 



Note that apparent frequency distribution of emergence does not change from the actual 
frequency distribution of emergence because number of cancellation with a decrease of 
in the system does not change. In the same manner, apparent probability distributions, 
apparent total occurrence rates, total change rates of flux amount, and apparent time 
scales are evaluated as 



QpAPP 

^f^(0) 

npAPP 
^ mrg 

dt 

apAPP 

dt 



(0) 



(0) 



npAPP 



2 / n{x)n{(p + x)m{4>,4' + x)dx 
n{x)l{(j), x)dx 





oo 



4't'ti 



k{x, (j) — x)dx 



2 / n{x)m{(j), x)dx. 

J<t>th 



(5.40) 
(5.41) 
(5.42) 
(5.43) 



dt 



= 2 



APP 
tot 



dt 



ON, 



APP 
tot 



dt 



spit 



ON, 



APP 



dt 



f 

0th 

poo 

J<l>f 



Ah 

oo 



POO 

n{(t)) I n{(f> + x)'m{(f>,(f> + x)dx 
Jo 

d(f) 



dcf) 



= 2 



n(0) / n{x)l{(f),x)dx 



n{(j)) I k{x, (f) — x)dx 

oo 



n{(j)) I n(x)m((f),x)dx 
<t>th 



dct) 



d(t). 



(5.44) 
(5.45) 
(5.46) 
(5.47) 



APP 
tot 



dt 



emrg 



n(0) / n{(f) + x)m{(l),(f) + x)dx 



(5.48) 
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APP 

tot 



dt 



APP 
tot 



dt 



mrg 



spit 



9$ 



APP 



dt 



oo 



4>n{(f)) I n{x)l{(j), x)dx 

Ah 



4>n{(j)) I k(x,(f) — x)dx 

POO 

4>n{(j)) I n(x)m((f), x)dx 

4>th L J 4>th 



(5.49) 
(5.50) 
(5.51) 



APP 



(f)n{(j))d(f) 



emrg 



/■oo 


poo 










4>n{(j)) I n{(j) 


+ x)m{(f), c 


p + x)dx 


d(j) 


•J (pth 


Jo 









T, 



APP 



(j)n{(f))d(j) 



4>th 



mrg 



4>n{(f)) I n{x)l{(j)^x)dx 
(j)n{(f))d(l) 



.APP 



th 
OO 



0th 



spit 



■4>~<t>t'ti 

4>n{(j)) I k(x, (j) — x)dx 



.APP 



(j)n{(f))d(j) 



4'th 



POO 


POO 








4>n{(j)) I n{x)m{q 


b, x)dx 


d(j) 











(5.52) 



(5.53) 



(5.54) 



(5.55) 



Next we evaluate total occurrence rate of apparent unipolar processes and total change 
rates of flux amount by them. There are four unipolar flux change events, namely unipolar 
appearance, unipolar increase, unipolar decrease, and unipolar disappearance. Figure [5?T] 
summarizes schematic pictures of them. 

Unipolar appearance corresponds to merging where flux contents of both merging 



patches are below (pth (Figure 5.1 1). Flux contents of both patches are below 0th- Sum of 
flux content is larger than (pth- Noting these conditions, total occurrence rate and total 



change rate of flux amount are deduced to 



dNi 



tot 



dt 



APP 



app 



20th r /•0th 

0th L^0-0th 



n{x)n{(j) — x)l{x, — x)d(p 



(5.56) 
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(a) Unipolar Appearance (b) Unipolar increase 




(c) Unipolar Decrease (d) Unipolar Disappearance 




Figure 5.1: Schematic pictures of unipolar flux change events, namely (a) Unipolar ap- 
pearance, (b) Unipolar increase, (c) Unipolar decrease, and (d) Unipolar disappearance. 
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dt 



APP 



app 



n{x)n{(f) — x)l{x, (j) — x)d(l) 



(5.57) 



Unipolar increase corresponds to merging where flux content of one merging patch is 



above (/)th but that of the other is below (pth (Figure 5.1 d). Total occurrence rate and total 



change rate of flux amount are deduced to in the same manner as 



dt 



APP 




r't>-4>t\i 










I n{x)n{q 


h — x)l{x, (, 


b — x)dx 


inc 




Jo 







d(j) 



+ 



(9$ 



tot 



dt 



APP 



2<Ath 



+ 



2<f>t 



n{x)n{(f) — x)l{x, (j) — x)dx 
n{x)n{(f) — x)l{x, (j) — x)dx 
n{x)n{(f) — x)l{x, cj) — x)dx 



(5.58) 



dcj). (5.59) 



Unipolar decrease corresponds to cancellation with a patch whose flux content is below 



dN,ot 


APP 


POO 




dt 


dec 







0th or splitting where flux content of one produced patch is below 0th (Figure 5.1 "). Total 
occurrence rate and total change rate of flux amount are deduced to in the same manner 
as 



n{x)n{(f) + x)m{x, (p + x)dx 

k{x,(j) — x)dx dcf) 

k{x,(j) — x)dx dip (5.60) 

n{x)n{(f) + x)m{x, + x)dx 
'4'th 

k{x, (j) — x)dx dip 
k{x,(j) — x)dx dcj). (5.61) 



/•20th 

+2 / n(0) 
>^0th 

+2 / ri(0) 



d^tot 


APP 


POO 


r<Pth 


dt 


dec 




Jo 



/'20th 

+2 / 0n(0) 

>^0th 
POO 

+2 / 0n(0) 



Unipolar disappearance corresponds to cancellation with a patch whose flux content 
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is below 0th or splitting where flux content of one produced patch is below 0th (Figure 



5.1 i). Total occurrence rate and total change rate of flux amount are deduced to in the 



same manner as 



tot 



dt 



tot 



dt 



APP 



mdisapp 



APP 



mdisapp 



2<At: 



.J(t>-<f>th 



n(0) 



n{x)n{4> + x)m{x, (p + x)dx 
k{x, (p — x)dx 

h 

n{x)n{(l) + x)m{x, + x)dx 

b 

k{x, — x)dx 



(5.62) 



(5.63) 



We evaluate lifetime of patches. Patch disappears in the analysis range when it cancels 
with a larger patch, merges with a larger patch or disappears by unipolar disappearance. 
Then disappearance probability, dP^f^p/dt{(f)), is deduced as 



(0) 



npAPP 
^ disapp 

dt 

1 apAPP f)pAPP npAPP 
"^^ tcnc I tmrg , mdisapp 

2 dt dt dt 
n(x)m{(f), x)dx + j n{x)l{(j)^x)dx + j k(x,(l) — x)dx 



(5.64) 



n{x)n{(f) + x)m{x, + x)dx (0th < < 20th) 

(20th < 0). 



n{x)m{(l), x)dx + I n{x)l{(f), x)dx 



The lifetime is an inverse of disappearance probability, namely 



r---(0) 

npAPP 
disapp 

dt 



(5.65) 



n -1 



(0) 
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n{x)m{(f), x)dx + j n{x)l{(j), x)dx + 
1 



+ 



n(0) 



n{x)n{(f) + x)m{x, + x)dx 



n{x)m{(j)^x)dx + / n(x)l{(f), x)dx 



4>th 

k{x, (p — x)dx 

(0th < < 20th) 
(20th < 0). 

There are two domains in the apparent lifetime. One is the range below 20th, where 
there is a unipolar disappearance by merging and cancellation and the apparent lifetime 
is different from the actual lifetime. Another is the range above 20th, where apparent 
lifetime is same as the actual lifetime. 



5.2 Discussions of Magnetic Patch Characters 
5.2.1 Flux Injection and Convective Effect 

Data set 1 and data set 2 have a different character in the variability of total flux amount. 
The total flux amount in data set 1 is almost constant during the observational period. On 
the contrary, that in data set 2 changes by ~ 50%. This difference is caused by large flux 
emergences, which may be injected from below the photosphere. Despite this difference, 
we find frequency distributions of flux content are the same power-law forms, which are 



also consistent with previous study by Parnell et al. (2009). It means that the injected 



flux is rapidly fragmented, namely surface processes maintain a frequency distribution of 
flux content. 
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5.2.2 Lifetime of the Patches 

We obtained the averaged lifetime of patches as 17.3 minutes in data set 1 and 23.2 minutes 
in data set 2. It is also found that lifetime has a power-law dependence as a number 



distribution on flux content with an index of -1.45 and -1.53, respectively (see Figure 3.7 



and Figure 3.10). The power-law distribution continues from ~ 70 minutes down to ~ 5 
minutes, which is near the time resolution of the data set, i.e. 3 minutes. The power-law 
distribution indicates that the averaged lifetime is not a typical one of the patches and 
that the value itself is not important. Further, we found that the corresponding power- 
law index of frequency distribution is ~ —2.45 and —2.53, respectively. The power-law 
index smaller than —2 means that patches with shorter time scale determine the obtained 
averaged lifetime. So the mechanism that determines the apparent lifetime is the process 
that is the most frequent one for their disappearance. From the results of Section 4.1, it 
is probably the unipolar decrease. See also Appendix B for the investigation of the flux 
amount of it. 



5.3 Discussions of Surface Processes of Magnetic Patches 

We did a statistical investigation of frequencies and change rates of flux amount by surface 
processes through an auto-detection code. The obtained results are as follows: 

1) Merging and splitting are much frequent than cancellation and emergence. 

2) Probability distribution of merging has only weak dependence on flux content, namely 
0.28 as a power-law index, and that of splitting has a similar tendency but rapid decrease 
around 0th- 

3) Frequency distribution of cancellation has a strong dependence on flux content, namely 
—2.70 as a power-law index. 

We speculate physical picture of flux maintenance in quiet regions from these results. 
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5.3.1 Comparison with Previous Studies of Bright Points 



Previous studies indicate that bright points in the photospheric and chromospheric lines 
(e.g. G-band, Call H, and Hq etc.) correspond to magnetic concentrations on the solar 



surface (Muller 1975 MuUer fc Keil 1983 Muller fc Roudier 1984 Berger fc Title 1996 



Ishikawa et al. 2007; Abramenko et al. , 2010). The time scales of merging and splitting of 



bright points are reported Berger & Title (1996); Berger et al. ( 1998b|[a ). They use high- 
resolution ground-based observation data taken by the Swedish Vacuum Solar Telescope 
on La Palma and investigate 534 bright points by an auto-tracking method. They report 
320 and 404 seconds as the mean occurrence times for merging and splitting of bright 



points (Berger et al. , 1998a), which are significantly shorter than our results, namely 



2000 seconds. We think that the difference of bright points may come from the fact that 



the bright point has a spatial scale of 200km (Utz et al. , 2009), which are shorter than 



our size threshold. The change of bright points occurs from the granular motion because 



the timescale of merging and splitting are the same order of that of granulation (Berger 



et al. , 1998a). On the contrary, our investigation may reflect the plasma motion larger 
than the granular motion because we set the size threshold of patches as 1.0 x 10^ km^, 
which is the typical size of granulation. The processes of bright points are controlled 
predominantly by the granular scale motion since the obtained time scales are similar 
(Berger et al. 1998a). On the other hand, those of our magnetic patches are controlled 
by the larger scale motion due to our setup in the size threshold. 



5.3.2 ProbabiUty Density Distributions of Merging and SpUtting 
in M-C equation 

We evaluate the probability density distributions of merging and splitting in a M-C equa- 
tion from the observational results in this section. Because the probability distributions, 
which we obtained in this study, are obtained by integrating them in the flux content 
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once, we have to put at least one assumption to evaluate them. 

The merging function l{x, y) is given as follows. We obtain the form as a probability 



distribution. By comparing Equation (5.41) and the result of Section 4.2.1 we obtain 



n(x)/(0, x)dx = po. 



mrg 



> 



(5.66) 



In the left-handside of this equation, the variable appears only in the /(0, y). We assume 
a simple form satisfying this relation namely. 



l{x,y) oc x^""^'^ 



(5.67) 



From the symmetry of l{x,y) = l{y,x), this relation deduces 



l{x,y) = Iq 



(5.68) 



Substituting it into Equation (5.66) and Equation (3.2) for n(0), we obtain 



X 



'dx 



(5.69) 



The upper value of integration is limited on the actual Sun and we put it as 0niax- We 
substitute the value obtained in the thesis , namely hq = 1.21 x 10^^^ Mx~^ cm~^. 



7 = 1.78, po,^ 



mrg 



2.56 X 10-4 s-\ l3-_ 



mrg 



0.28, 



10^^ Mx, 0th = 10^^-^ Mx, and 



10^^ Mx, and obtain Iq = 7.24 x lO^^ cm^ s'^. 



Next we evaluate the probability density function of splitting. The apparent frequency 
distribution of splitting has more complicated form than that of merging because splitting 
into patches below the detection threshold changes actual frequency distribution more 
than merging. It is difficult to evaluate the actual frequency distribution from the apparent 
one. The strategy we choose here is that we assume the form of splitting function i.e.. 
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k{x, y), and evaluate the apparent frequency distribution from that. Then we can see the 
difference between the evaluation and the observational result. From our observations, the 
probability distribution of splitting events dP^^^{(j))/dt is suggested to be independent 
of the parent patch flux: 

— — — (0) = ko = constant (for all 0). (5.70) 



This claim is observationally supported at least in the range > 0th (see Figure 4.4). 
The drop off below 0th is discussed immediately below. If the splitting ratio between the 
daughter patches is randomly determined i.e., 

d 

— [fc(x,0-x)] = (for < x < 0), (5.71) 



then we obtain, 



k{x,y) = ^. (5.72) 
X + y 



When the flux content of the daughter patch is below 0th, such events are not recognized 
as a splitting event in our procedure. The probability distribution will be given as 



^pAPP /■<f>-<t>th 2rh.u 

^P'*.(0)=/ k{x,(l)-x)dx = ko{l- (5.73) 



dt 



The black and blue dashed curves in Figure [4^ indicate analytical curves with ko = 1.0 x 
10^^ s^^ and /cq = 5.0 x 10^^ s^^, respectively. This curve fits the drop of the observational 
line well, which supports the above assumptions. We obtain k{x,y) = ko/{x + y) as a 
splitting function, where ko = 5.0 x 10^^ s"^. 



5.3.3 The time Independent Solution of Splitting Process 

Since our observations show that merging and splitting are much more frequent than 
emergence and cancellation, it is suggested that the former two have much influence of 
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the maintenance on the power-law distribution. Though we have not find the solution of 
M-C equation with both merging and splitting, we show the time-independent solution 
only with splitting. Numerical calculation of M-C equation may be useful to find the 
solution with splitting and merging in the future works. 

The frequency of emergence, merging, splitting, and cancellation are represented by 
S{(j)), l{x,y), k{x,y), and m{x,y) respectively. We use the equation only including the 
splitting terms by setting 5(0) = 0, l{x, y) = 0, m{x, y) = 0, namely 



n(x)k{(l), X — (j)) dx — / n{(j))k{x, (j) — x) dx. (5-74) 







dn{(f)) 
dt ^ 

Here the splitting process, k{x,y) is only taken into account. We obtained the function 



form of k{x,y) in Section 5.3.2 and obtained k{x,y) = k^/^x + y). After substituting it 



into Equation (5.74) and differentiating with 0, we obtain 

d^n{(j)) ko d 



d(t)dt 



h2 



[<PM<I>)] ■ (5.75) 



This equation has a time-independent solution n{(j)) oc 0"^. This power-law index of the 
flux content is in good agreement with the observational result, 7 = —1.8. It is possible 
that n(0) has a solution with a power-law function by using both the splitting and merging 
processes in the above analysis though it is not simple task due to its mathermatical 
difficulty. One further investigation is to investigate time evolution of the frequency 
distribution of the flux content during decay of active regions and frequencies of merging 
and splitting. If the merging and splitting make the power-law distribution of the flux 
content, we will observe the change of the power-law index to —1.8. 

This scale-free distribution comes from constant dependencies of splitting, namely 
constancies on parent flux content and daughter flux content. These constancies may 
come from splitting dominated by convection or flux tube instability. In order to distinct 
these models, we propose two future observations: One is an investigation of convection 
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flow field by using Hinode/FG data. The other is that of fine-scale structure in the flux 
tube and its dynamics by using Hinode/ SP. 



5.3.4 Interpretation of Frequency Distribution of Cancellation 
and Emergence 

We discuss why the slope of frequency distribution of cancellation is so steep compared to 
those of merging and splitting in this section. We suggest a relationship model among the 
frequency distributions of the flux content, cancellation, and emergence. The constant 



velocity and randomness of direction of motion are assumed here. Figure |5.2| shows a 
schematic view of this model. We put a power-law distribution of flux content as 



n(0) = no 



(5.76) 



The power-law index, 7, is derived as 1.5 < 7 < 2 by our observation and the previous 



studies (Parnell et al. , 2009, Zhang et al. , 2010). The maximum flux content in the system 
(0max) is assumed to be much larger than the minimum (0min) in the following discussion. 
We calculate N{(f)), a total patch number with a flux content larger than (p, by integrating 



from (h to 



as 



iV(0) 



n{(j)) d(j) 



no<f)o f 4> 
7-1 



-7+1 



(5.77) 



Schrijver et al. (1997) evaluated a collision rate of opposite patches from a total patch 



number with assumptions of a constant velocity and a randomness of patch motion along 
a network. They obtained the collision frequency, u, as 



4VP 



(5.78) 



where Vq, p, and A^^ mean a typical velocity of patches, a number density of network 
cell, and a number density of patches, respectively. We multiply 1/2 taking the double 



CHAPTER 5. DISCUSSIONS AND INTERPRETATIONS 



88 



counting into account. We obtained that the frequency of merging and splitting are much 
larger than that of cancellation. It can be deduced that the frequency distribution of 
flux content is maintained rapidly by merging and splitting compared to the time scale 
of cancellation. This enables us to treat that the number density of patches is time 
independent in the evaluation of cancellation and apply the same analogy to the number 
density expanded in the dimension of the flux content. The collision frequency of opposite 



polarities, dN{(f))/dt 



, is evaluated as 

col 



dN{(f)) 



dt 



4(7-l)2v^ 

col 



(5.79) 



Note that this total collision number becomes time-independent with an assumption of 
maintenance of a power-law distribution of flux content. We assume that the total number 
of cancellation, dNcnc{4>) 1 9t, equals the total number of collision events of the opposite 
polarities, dNcoi{4>)/dt, namely, 

dNU<t>) giVcoi(0) 

The frequency distribution of cancellation, dncnc{<P) / dt is given by differentiating equation 
with 0, namely 



dt d(pdt 2(7-l)yp 

This assumption means that there are no scripting patches in the collision of opposite 
polarities once they collide, which is difficult to check and we justify this assumption from 
the comparison of the obtained frequency distribution of cancellation in this model and 
that in the observation. Based on the above discussion our estimation predicts a frequency 
distribution of cancellation as dricnc/dt ~ -9.4 x IQ-^^ ^ (0/i.o x 10^^ Mx)"^'^'^^"'^^ Mx'^ 



cm ^ s ^ where we adopt p = 1.0 x 10 cm ^ (Hagenaar et al. , 1997), Vo = 1.2 km s ^ 
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(Section 3.4), no = 6.8 x 10^'^^ Mx^^ cm^^ (Section 3.3), and 7 = 1.93 from data set 
2. The observed frequency distribution of cancellation for data set 2 is ducnc/dt = 
— 1.29 X 10""^^ X (0/1.0 X 10^^Mx)~^ ''^^°'^^. The absolute values of the distributions are 
in consistent within a factor of 2 and the power-law index is in consistent within the 
l(T-error. The difference of absolute value may be explained by the fact that the patches 
have the finite size and it increase the rate of the collisions. 

We also evaluate a frequency distribution of emergence with assumptions of small 
amount of flux supply from the outside of the system and re-emergences of canceled fluxes. 
These assumptions lead to the relationship in which a flux distribution of emergence nearly 
equates that of cancellation, 

9raemrg(0) _ <9ncnc(0)_ ^^O^o'/'O Z' ^ 



dt dt 2(7-l)^V</'o 



Thornton & Parnell (2011) reports a power-law frequency of emergence as c^'^emrg/ ^ 
1.5 X 10"''° X (0/1.0 X 10^® Mx)"^''' Mx"-' cm"^ s"-'. The steepness of the estimated power- 
law distribution, -2.48±0.24, is in good agreement with the observational results, -2.7. It 
supports the idea of recycling of flux content in a quiet region. The discrepancy in absolute 



value between our estimation and Thornton & Parnell (2011 ) may be due to the difference 



in patch detection methods. It is reported (DeForest et al. , 2007) that there is a two- 



orders-of-magnitude enhancement in the number of detected patches by the clumping 



method (ours) over the downhill one (Thornton & Parnell, 2011). To investigate this 



difference, we will investigate the frequency distributions of emergence and cancellation 
in the same data set with the same detection method as a future work. As for the direct 
investigation, the flux recycling is difficult to be determined because we cannot observe 
the magnetic structure below the solar surface. One possible investigation is to make pairs 
of emergence and cancellation by the correlation time and the distance between them. 

Further we calculate the time scales of cancellation and emergence to see their de- 
pendences on the detection limit, which is reported from the observation as that they 
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become much shorter with the higher resolution. The total flux amount in the system is 



calculated from Equation (5.76) as 



tot 



4>n{(j)) d(j) 



y^o0o f 4>: 
2-7 



-7+2 



(5.83) 



Since 7 < 2, this result shows that the total flux is dominated by patches with larger 
flux content. On the other hand, we obtain a total flux loss amount by cancellation, 



9$tot(0)/<9t|(,nc5 froi^ Equation (5.81) as 



tot 



dt 



dn{(t)) 



dt 



v^nlcpl 



2(7-l)(27-3)^ 



-27+3 



(5.84) 



The total flux supply by recycled emergence is also evaluated in the same manner from 



Equation (5.82) as 



9$t 



dt 



emrg 



dn{(f)) 



dt 



Vonl(f)l 



emrg 



2(7-1)(27-3)Vp V 00 



Then time scales of cancellation and emergence are evaluated as 



-27+3 



(5.85) 



emrg 



$tot/(9$tot/5t|,,J = $tot/(9<l>tot/9t| 

2(7-l)(27-3) 



(2 - -f)vono(j)oy/p V 00 



-7+2 



emrg/ 
27-3 



(5.86) 



The actual value of 7 was obtained as 1.8 — 1.9 in this study, which means that 27 — 
3 = 0.6 — 0.8. The positive dependence on the detection limit, which is represented by 
0min here, means that the flux replacement time scales becomes shorter with the smaller 
detection limit. This result is qualitatively consistent with the observational fact that the 
flux replacement time scales by emergence and cancellation becomes shorter with higher 



resolution in the previous studies (Martin et al. , 1985; Schrijver et al. , 1998, Hagenaar 



2001). 
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Figure 5.2: Relationship model among the frequency distributions of the flux content, 
cancellation, and emergence with 1.5 < 7 < 2. The solid arrows indicate mathematical 
relations, while the dashed arrows indicate relationships with some assumptions. 
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5.4 The Interpretation of Flux Maintenance in Quiet 
Regions 



We summarize our interpretation from the discussion here. Figure [573] shows a schematic 
picture of the interpretation. First, a power-law frequency distribution of flux content 
is rapidly maintained by merging and splitting, which is supported by the result of a 
comparison of change rate of flux amount by each activity. In addition to this, we have 
found that a power-law frequency distribution with an index of —2 is a time-independent 
solution of splitting. Cancellation is caused by convective motion and frequency distribu- 
tion of cancellation should naturally become a steep power-law distribution. Emergence 
is interpreted as a re-emergence of submerged flux by cancellation, which is consistent 
with a steep frequency distribution of cancellation in our study and emergence reported 



by Thornton & Parnell (2011). One of the possibilities but the important point of this 
model is that the apparent flux transport by emergence and cancellation becomes much 
more drastic with higher resolution even when there is no transport. 
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Flux Injection from below 
photosphere and/or AR 



Rapid maintenance by 
splitting and merging 



Dist. of Flux Content 

Index: 
-1.8-1.9 



Dist. of Emergence 



Index: -2.7 



Patch motion by 
convection 



Dist. of Cancellation 

Index: 
-2.4-2.7 



Re-emergence of 
canceled flux 



Figure 5.3: Schematic picture of our speculation of flux maintenance in quiet regions. 
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Chapter 6 



Conclusion and Future Works 



6.1 Conclusion 

Our conclusion is that frequency distribution of magnetic flux in quiet regions is main- 
tained by surface magnetic processes. The frequency distribution of the flux content has 
a powe-law dependence with an idex of —1.8 — 1.9 which continues to 10^"^'^ Mx at least. 
The patches with this distribution move randomly and with a typical velocity, which may 
be caused by convective motion on the solar surface. The frequency distribution of can- 
cellation has a steep power-law distribution with an index of —2.48 ± 0.24. This steep 
distribution is explained by a random convective motion and collisions of patches. We 
speculate that the patches submerge through cancellation and re-emerge to solar surface 
as an emergence, which is supported by comparing the obtained frequency of cancellation 



and that of emergence reported by Thornton & Parnell (2011). Splitting and merging are 



caused with an timescale of ~ 33 minutes on the solar surface. These timescales have 
only a weak dependence on flux content, —0.28 as a power-law index for merging and 
no significant dependence for splitting. These weak dependence may cause a power-law 
distribution with an index of ~ —2. 

In this thesis, we develop an auto-detection code of surface processes and apply it to 
two different data sets of magnetograms obtained by Hinode spacecraft. In Chapter 4 
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of our thesis, we investigate relationships between flux content and size, proper velocity, 
and lifetime of patches statistically. The proper velocity is found to have only weak 
dependence, —0.28 as a power-law index. We treat the proper velocity of patches as 
constant in the later discussion from this result. It is found that small patches determine 
the averaged life. It suggests a possibility that the apparent fragmentation processes to 
undetected scales, namely cancellation and splitting involving undetected patches, are 
much frequent than actual disappearances. On the other hand, there is an upper limit of 
lifetime, ~ 60 minutes. These results indicates the dissociative nature of patches in the 
analyzed flux range. 

The measurements of probability distributions of surface processes are done in Chapter 
5. It is found that splitting and merging are much more frequenct than cancellation and 
emergence. We further investigate frequencies of merging, splitting, and cancellation. 
Probability distribution of merging has only weak dependence on flux content with an 
index of —0.28. That of splitting is also almost constant but there is rapid decrease 
around the detection limit. Frequency distribution of cancellation has a steep power- 
law distribution with an index of —2.48 ± 0.24. The obtained results are converted to 
probability density distributions. Analysis of M-C equation shows that splitting has a 
time-independent solution of power-law distribution with an index of —2. On the other 
hand, the steep power-law of frequency distribution of cancellation can be explained with 
a constant proper velocity evaluation, which we interpret they are advected by convective 
motion. 
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6.2 Future Works 

6.2.1 Extending Studies of Magnetic Structures on the Solar Sur- 
face 

In this thesis, we estimate each term in M-C equation, i.e. l{x,y), k{x,y), and m{x,y), 
in an indirect manner under an assumption of the symmetry in one of two dependent 
variables (See Section 5.3.2 for detail). The main reason is a lack in numbers of events to 
obtain enough statistics. In order to do a direct measurement, two orders of magnitude 
larger number of events are necessary. At this moment, such data is unavailable in the 
Hinode database. The full Sun observation data of magnetic field by SDO is what we 
should work out in the near future. 

The discussion in this thesis is based on the M-C equation that only has the time 
and magnetic flux as the independent variables, while the actual processes on the solar 
surface are much more complex with many degrees of freedom. For example, the filling- 
factor of the magnetic field is an important factor that controls the time scale and the 
rate of flux re-distribution in a splitting event. If the filling-factor is smaller in larger 
patches, then the convective flow can split them more easily. This may explain our result 



of non-dependence of splitting frequency on the flux content. Convective collapse (Parker 



1978 , Nagata et al. , 2008 , Fischer et al. , 2009 ) in larger patches may explain this relative 



difference in the filling factor. Using the Hinode/ SP observations, the detailed information 
in each splitting event would become available and will make it possible to evaluate this 
model. 

The collaborative investigation with numerical calculations is useful. We suggest the 
qualitative picture in this study but we need to confirm whether it is possible in quanti- 
tative sense. It is achieved by solving the M-C equation numerically by using our mea- 
surements. 



Meyer et al. (2011) report their two-dimensional numerical simulations based 



on the modeling (Parnell, 2002). They investigate equilibrium state by injecting fluxes 
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and controlling the splitting frequency as a free parameter. A power-law distribution of 
flux content is maintained with an input of a steep power-law frequency of emergence. 
The induced frequency of cancellation becomes a steep power-law one, which is consistent 
with our model. 

Recent progresses of numerical simulation enable more realistic simulations on the 
solar surface, namely radiative magnetohydrodynamics simulation. Although it is still 
difficult to set the numerical domain large enough for a statistical investigation, such 
calculation would be very helpful to determine the physical interpretation of frequencies 
of surface processes of magnetic patches. 



6.2.2 Contributions to Studies of Other Solar Phenomena 



We think there are two solar phenomena that are related to the results of this thesis, 
namely the surface flux transport and the nature of X-ray bright point. 

The flux transport to the solar poles is an essential ingredient of the fiux-transport 



dynamo model (Wang et al. , 1991). The major part of the flux in this model is provided 



by the diffusion of active regions. So it is necessary to extend to cover the active regions 
beyond our current analysis limitted to quiet regions. It may also be necessary to cover 
much smaller scale including the intranetwork field. 



Transient Horizontal Magnetic Field (THMF) (Figure 6.1 Ishikawa & Tsuneta, 2009 



Ishikawa et al. , 2010). is the horizontal field with a short lifetime of 1 — 10 minutes 
and a flux content smaller than lO^'''^ Mx. Ishikawa & Tsuneta (2010) reports a three- 



dimensional view of THMF, that THMF is an emerging flux tube within the internetwork. 
THMF can be the recycled magnetic flux, which is supported by the fact that the flux 



amount of THMF, 1.8 x 10^^ Mx day ^ for the entire surface (Ishikawa & Tsuneta, 2010) 



is as large as the small emerging bi-pole detected in (Thornton & Parnell, 2011), 3 x 10 



125 



Mx day ^ for the entire surface. Of course we need to investigate the spatial and temporal 
relationship between cancellations and THMF as the future works to give a conclusion to 
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this phenomena. 

For the future investigation of these small field in the internetwork, the next high- 
resolution spacecraft mission, Solar- C, will be a great help. It is planned to have a higher 
resolution by one-order of magnitude in magnetic flux content. We may be able to set the 
detection limit as 10^^'^ Mx, which corresponds to one-order of magnitude in internetwork 
field domain. It will be possible to do the statistical survey of such internetwork field based 
on the auto-detection of patches in Solar- C era. 

The other extension of our analysis is to investigate the latitudinal difference of oc- 
currence of cancellations. We need to have the patches submerge or at least disappear at 
the high latitude for the polar reversal in solar cycles, both of which can be done only 
through cancellations. So there is a possibility that there is a latitude dependence of 
cancellations. Recent paper, Shiota et al. (accepted), reports another interesting result 
that the decrease in the net flux is caused by a decrease in the number and size of the 
large flux concentrations (> 10^® Mx) as well as the appearance of patches with opposite 
polarity at lower latitudes. Our method of analysis will be also helpful to decide whether 
these large patches are convected as a small patches and maintained by mergings in the 
polar regions or convected as the large patches through the quiet regions. 

The energetic events on the solar surface are related to magnetic activities as mentioned 
in Chapter 1. X-ray Bright Points (XBPs) are the brightenings seen in the X-ray images 



Golub et al. (1976). One model of XBPs is a magnetic reconnection in the higher solar 



atmosphere between different flux systems (Figure 1.4). This model is suggested from 



the fact that most XBPs correspond to cancellations (80%; Harvey, 1997), which are 



thought as submerging loops after reconnection (lida et al. , 2010). We obtain an energy 



distribution of BPs from Equation (5.81) and these relationships as 



-27 — anp+2 

,2^2 / Z? \ agp 



^"BP(^) = ^"^o'^o (^) " . (6.1) 



dt (7 - l)aBPy/pEo \Eo^ 

We set two assumptions here. One is that all cancellations coincide magnetic reconnection 
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Figure 6.1: A schematic picture of network and internetwork field in quiet regions. The 
blue (red) color in background shows upflow (downflow) motion in supergranule time 
scale. The green arrows indicate the horizontal convective motion. Internetwork field 
is located in the supergranular cell. It is thought to be disturbed by smaller convective 
motion, namely granular cell and turbulence in the cell. The network field is a magnetic 
field swept to the edge of supergranular cells. The one at the conjunction point has a 
larger size and is more stable. 
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in higher atmosphere, which is supported by the observational results in the Yohkoh era. 
The other is that there is one-to-one relationship of a power-law distribution between 
canceling flux content and released energy, namely dn^p/dt = Eq {(p / (po)"'^^ . Rough 
estimation, 7 = 1.8 and Obp = 1.5, deduces a power-law index of BP energy as —2.1. 
The power-law index steeper than -2 means that the smaller XBPs has a significant 
effect for the total energy release. We want to note that the scale-free nature of energy 
distribution of BP comes from scaleless of flux content in this model. There are some 



papers reporting the power-law index of the energy release of XBPs Crosby et al. (1993); 



Shimizu & Tsuneta (1997); Krucker & Benz (1998); Parnell & Jupp (2000); Benz & 



Krucker (2002); Aschwanden & Parnell (2002). Some show that steeper than -2 but most 



show the power-law of —1.5 1.7, which are flatter than our estimation. However, it 

should be noted that their target is the XBPs in the loop structure, namely not canceling 
XBPs. We speculate that the appropriate model for the XBPs in the loop structure 



is the entangled flux tube suggested by Parker (1972), where the energy is released by 
the component magnetic reconnections of nearly paralleled field. Thus we expect the 
different distribution of released energy as XBPs from our estimation. The investigation 
distinguishing canceling XBPs and XBPs in the loop structure as the future works. 
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Appendix A 

Flux Imbalance in Merging and 
Splitting: Does It Affect Frequency 
Distribution of Cancellation? 



Because merging and splitting are much more frequent than cancellation in our result, 
there is a possibility that flux imbalances during merging and splitting have significant 
effects in the evaluation of the frequency distributions of cancellation. We examine this 
possibility by comparing the frequency distribution of flux imbalances during these events. 

First we investigate flux imbalances during merging and splitting in data set 1. Figure 
A.l and Figure A. 2 shows scatter plots of flux content in parent patch and daughter 



patches. We see that flux amount seems to be conserved well in the both processes from 
these figures. 

Next we investigate frequency distributions of flux imbalance of these processes. The 
flux imbalance is defined as difference between sum of flux amount in parent patches and 
that in daughter patches, namely 

"^parent ^daughter 

i=0 j=0 
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Figure A.3 shows the result. The solid histograms indicate that for merging and the 
dashed ones for splitting. A green vertical line shows the detection limit of this study. 
First, we see that flux imbalance is less than the detection limit. It shows well-balanced 
flux content through merging and splitting. Second, magnitude of frequency distribution 
in the range above detection limit is larger than that of cancellation, which is shown by 
black histogram, by one — two orders of magnitude in the analysis range. We therefore 
conclude that flux imbalances in merging and splitting do not have a significant influence 
on frequency distribution of flux content. 




Figure A.l: A scatter plot of the flux content before and after merging processes in data 
set 1. 



APPENDIX A. FLUX IMBALANCE IN MERGING AND SPLITTINQOG 




Figure A. 2: A scatter plot of the flux content before and after splitting processes in data 
set 1. 
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Frequency Distributions of Flux Difference 
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Figure A. 3: Frequency distribution of flux imbalance in merging and splitting processes. 
The red and blue solid histograms indicate that for merging of positive and negative 
polarity and the dashed ones for splitting. The black solid histogram shows the frequency 
distribution of cancellation in data set 2. 
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Appendix B 



Unipolar Processes Observed in This 
Study 



We found many unipolar events, namely unipolar increase, appearance, decrease, and 
disappearance in our data sets as shown in Table 4.1. Here we show examples of unipolar 
processes observed in this study to support the assumption that these events are interac- 
tions with patches below the detection limit, which is put throughout the thesis. 

Although it is often mentioned that plenty of unipolar events are observed in magne- 



tograms, there are few studies about these unipolar events. Lamb et al. (2010) shows one 
example of unipolar appearances as a coalescence of patches below a resolution limit by 
using magnetograms with different spatial resolutions. 

We interpret the unipolar events as interaction with patches below the detection limit. 



Figure B.l shows examples of unipolar processes detected in data set 1. The background 
shows the magnetic flux density and the red (blue) contours indicate the positive (nega- 
tive) patches with the detection condition. Figure B.l[ a) shows a unipolar increase. There 



is a negative patch at 1:46UT indicated by the blue contour with the red arrow (patch 
A). There are two small patches without contours (patches B & C) in the east-side of the 
detected patch. They are gradually absorbed to the patch A and they are recognized as 
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one large patch at 2:01UT. The patch stays stable and recognized as one magnetic massif 
after the coalescence. Because the patches B & C are below our detection limit, this event 
is recognized simply as an increase of the flux content in the patch A, i.e., a unipolar in- 
crease event. Figure B.l[ b) shows a unipolar appearance of positive patch. In this case. 



there are three small patches below the detection limit at 2:29-2:31UT. They converge and 
coalesce to one large patch at 2:33UT. It is difficult to see the inner structure of the patch 
even with the high spatial resolution of Hinode spacecraft. Figure B.l[ c) shows a unipolar 



decrease. One detected positive patch is near the center of view at 1:41UT. The east part 
of the patch changes its shape and one patchy structure can be seen there. It gradually 
grows during 1:42— 1:45UT and is finally separated at 1:49UT. The separated patch is 
small and it is below the detection limit. Figure B.l[ d) shows a unipolar disappearance 



of the negative patch. One negative patch is detected at 2:19UT. However, a dipped 
shape is suddenly formed in the outline of the patch at 2:27UT. The patch is divided 
into two small negative patches below the detection limit at 2:29UT. These four events 
are explained as interactions with patches below the detection limit, which is the same 



result as (Lamb et al. , 2010). Through our interpretation, the unipolar events should be 
counted for splitting events. This leads to the underestimation of the number of splitting 
events and is the cause of an apparent drop in the probability distribution (Figure 4.4). 
Such infiuence of the detection limit is taken into account in our analysis (Section 5.1.3). 

Further, we investigate the fiux amount changed by the unipolar processes in data set 
1. It should concentrate around the detection limit based on our interpretation. Figure 



B.2 shows the histogram of the total fiux amount changed by unipolar processes. We 
can see that all four histograms concentrate around the detection limit, namely the green 
vertical line. This result supports the interpretation that the unipolar processes are the 
interaction with the patches below the detection limit. Because we take the infiuence of 
the detection limit, the power-law index induced in the thesis may not be affected by 
these unipolar events. 
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Figure B.l: Examples of apparent unipolar events, namely, (a) unipolar increase, (b) 
unipolar appearance, (c) unipolar decrease, and (d) unipolar disappearance. The back- 
ground shows the magnetic flux density obtained by Hinode/NFl. The red (blue) con- 
tours indicate positive (negative) patches detected with our threshold. The field of view 
is 9.6" X 19.2" for all images. 
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The further quantitative investigation of unipolar processes is needed to give a con- 
clusion for the discussion here. We believe that the auto-detection with the different 
conditions will be useful at this point and it should be done in the future works. 
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Flux Content of Unipolar Processes 
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Figure B.2: Histogram of flux amount changed by unipolar processes in data set 1. The 
solid lines indicate increase processes and the dotted lines indicate decrease processes. 
Colors indicate the polarity, namely red for positive and blue for negative. 
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